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Die Erwärmung in metallischen Kreislochscheiben bei 
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radialem Heizstrom und temperaturabhängigem Leiterwiderstand 


Von H. BucHHorz - ° 


Mittels der Laplace-Transformation wird in der vorliegenden Arbeit im wesentlichen der Er- 
wärmungsvorgang in einem Hohlzylinder berechnet, wenn u.a. die Erwärmung des Zylinders durch einen 
‚in radialer Richtung fließenden elektrischen Strom erfolgt. Im Hinblick auf die kurze axiale Länge des 

Zylinders ist nur die Abhängigkeit der Temperatur von der radialen Koordinate o von Belang. Da die Tempe- 
 raturabhängigkeit des spezifischen elektrischen Widerstandes berücksichtigt wird, wird die transformierte 
Differentialgleichung des Erwärmungsvorganges inhomogen von der zweiten Ordnung, so daß in der Lösung 
neben den Zylinderfunktionen auch die Lommelsche Funktion auftritt. Der Parameter der Zylinder- 
funktionen wird rein imaginär, aber nur klein. Für die Nullstellen der Koeffizienten-Determinante und für Be 
die Zylinderfunktionen und die Lommelsche Funktion bei rein imaginärem Parameter werden Näherungs- 
ausdrücke angegeben. f ; 


. Using Laplace transformation, the paper calculates the heating in a hollow cylinder e. g. due to an electric 
current in radial direction. T'he axial length being small, it is only the dependence of temperature on the 
radial coordinate that matters. Since temperature dependence of electric resistivity is taken into account 
the transformed differential equation describing the heating process is inhomogeneous of second order; the 

Br .. solution, therefore, contains the Lommel function besides the Bessel functions. The parameter of the Bessel 

2 7 functions will be purely imaginary but small. Approwimate expressions are found for the zeros of the deter- 
minant of the coefficients, as well as for the Bessel functions and the Lommel function for purely imaginary 
values of the parameter. 


von 


B HacTosmei pa6oTe C MOMOIIBIO Ipeo6pasoBaHunn JlamMlaca BbIUHCJAeTCcH Ipolecc Ha- 
TpeBaHuA B IHOJOM UNJIHHAPe, HAIPHMep, KOTNa HATPeBaHHe INJIMHAPA IIPOHCXOAHT UON Nei- 
CTBHEM IANEKTPHYECKOTO TOKA, KOTOPbIH IPOTeRAaeT B PANHAAIbHOM HanpaBıteHun. BBuANy Toro, 
YTO BbICOTA UNMJIHHAPA OYeHb Mala, HMeeT 3HAYeHHE TOJIBKO BABHUCHMOCTB TEMIUEPATYP&I OT pa- 
AHAAbBHOÜ KOOPAHHATE 0. Tak KAaRK YUHTBIBACeTCA 3ABHCHMOCTb YAEIBHOTO IANEKTPHYECKOTO 
COUPOTUBJIEHNHA OT TEMIIEPATypbI, TO IIpeo6pasoBaHHoe AuhdhepeHmmanbHoe YpaBHeHue IIPO- 
ıecca HaTpeBaHun ÖYyNEeT HeONHOPONHEBIM BTOPOTO Hopamka. B pemeHHe KpoMe IHJIHHAPHYE- 
CKUX PyHKUHÜ BXONUT Tarıke hPyHruna JIomMesn. llapamerp unsmHapmyecKof PYHRIMHM 
ÖyAHeT YMCTO MHUMbIM U MAJIEHLRUM. Jia KOpHeä ofpenesnTenan KoaßPHLUMeHTOB, MIA IM- 
AUHAPHUYecKuX PyHRUna u num Pyukumnm JIoMMe,ıa Opu YMCTO MHMMOM IapaMerpe yRa3aaHbl 
IPHÖJIBKEHHA. 


N \ u 2% a ll nn Sk 


1. Die besonderen Annahmen der Aufgabe 


Az 


Den einzelnen Elementen in den thermoelektrischen 
Batterien wird neuerdings [1], [2] häufig die Form dünner, 
plattenähnlicher Kreislochscheiben gegeben, die bei ko- 
axialer Anordnung gemäß Bild nur durch schmale Luft- 
zwischenräume voneinander getrennt sind. Die inneren 
und äußeren Ränder der Scheiben sind jeweils oben oder 
unten zusammengefaßt, und durch besondere dort einge- 
baute Heiz- oder Kühlplatten läßt es sich mit ausreichen- 
der Annäherung erreichen, daß der Außenrand der Scheibe 
mit dem Radius o=b vom Augenblick der Inbetrieb- 
nahme an, d.h. für alle t > 0 die gleichbleibende Tempe- 
ratur 9, und der innere Rand mit dem Radius o = a die 
feste Temperatur 9, behält. 

Wird das Zylinderkoordinatensystem zu der Scheibe 
nach Art des Bildes angeordnet, so läßt sich die Tempera- 
turverteilung in den einzelnen Scheiben bei deren Vielzahl 
und wegen der geringen Dicke 2 d der einzelnen Platten in 
ausreichender Annäherung als nur von o und f abhängig an- 
sehen und also etwa durch die Funktion (eo, t) beschreiben. 
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Es sei dann zum Zeitpunkt t = 0 
%K,0)=0 für asosb ......... !rr0E 


Außerdem ist wegen der Kühl- und Heizplatten längs des Umfanges der Scheibenränder 0 = a 
unde=b | 
Mah=d, » =: Dr ee 
HD =d 5, He ren 


mit a<b und 9,< d,. Wird dann zum Zeitpunkt {=0 auch der Stromkreis des Thermo- 
elementes geschlossen, so entwickelt sich zunächst einmal in den Kreislochscheiben ein Wärme- 
strom infolge dieses gleichbleibenden Temperaturunterschiedes zwischen den heißen und kalten 
Scheibenrändern. Der in den Scheiben in radialer Richtung fließende Strom I im Stromkreis 
verursacht jedoch nach dem Gesetz von JouULE einen weiteren, zusätzlichen Wärmestrom. 

Um den unter diesen Voraussetzungen zum Zeitpunkt f= 0 einsetzenden Erwärmungs- 
vorgang analytisch zu behandeln, führen wir die folgenden Konstanten des Wärmestromkreises 
ein: Es werden bezeichnet mit 

[6 die spezifische Wärme des Scheibenmaterials mit der Dimension [ce] = Ws/(cm? - grd) 

A die Wärmeleitfähigkeit mit der Dimension [A] = W/(cm - grd). 


Ferner sollen bedeuten: 

2d die Dicke der einzelnen Scheiben mit [d] = em, . 

0, den spezifischen, elektrischen Widerstand bei der Ausgangstemperatur = 0 mit 
[o,] = Ohm - cm, 

& den Temperaturkoeffizienten des Scheibenmaterials mit [a] = 1/grd (Zahlenwert: 
a» 4.10), 

I den gesamten, gleichmäßig über den Leiterquerschnitt verteilten elektrischen Strom 
in jeder einzelnen Scheibe in A und also 


I/(2ro.2d) die Stromdichte in der Scheibe in A/cm? im Abstand von der Achse der Scheibe 
füreinomitasosb. 


Von den drei Materialkonstanten c, A undo=o,(l1 + ad) soll also im Rahmen dieser 
Arbeit nur der spezifische Widerstand o als temperaturabhängig angesehen werden. Dies läßt 
sich im Hinblick auf die Kopplung zwischen o und A über das WIEDEMANN-FRANZsche Gesetz 
nur deshalb vertreten, weil die Änderungen der Wärmeleitfähigkeit gegenüber denen des spezifi- 
schen Widerstandes mit der Temperatur lange nicht so stark ins Gewicht fallen. Streng genom- 
men bringt dann diese Zunahme von o mit der Temperatur bei konstanten Werten von 9, — Ö, 
allerdings auch wieder eine Anderung von I mit sich. Es muß also hier noch zusätzlich voraus- 
gesetzt werden, daß die übrigen Stromkreiswiderstände, die ja wegen der weit geringeren Höhe 
der Temperatur kaum in ihrer Größe beeinflußt werden, die Zunahme des Eigenwiderstandes 
der Scheiben nahezu völlig überdecken. 

: Dann gilt bekanntlich nach dem Energieprinzip für (o, t) die folgende Differentialglei- 
chung 


e ae BURN ae 2m 2 =9-(1+0:%6,)- B Errer)] Me: 

mit den weiterhin immer wiederkehrenden zusammengesetzten Abkürzungen: 
9,=ojA-(Iärnd® mit [9,]=grd und [a-@,l=1...... 0a) 
=eis=i(609,)% mt plelkls1i ud ief 12.202260 
Ale=yR mit [y2] = CB ne 


Damit kann die Gl. (2) auch in der Form geschrieben werden 


,„ Ro, 1 
EN 


‚ao © 9%Xo,t) 
00°? do y? ot 


. + & 0, . %o, I) == 


Ge 


Um diese inhomogene, partielle Differential 
Anfangs- und R 
werden. 


\ 5 N gleichung der Wärmeleitung unter den angegebenen 
andbedingungen aufzulösen, soll hier die LAPLACE-Transformation herangezogen 


Mi 
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2. Die Auflösung der inhomogenen transformierten Differentialgleichung 


Wie es vielfach wegen der einfachen Schreibweise üblich geworden ist, möge in dieser Arbeit 
die LAPLACE-Transformierte von %e, f) durch [8, S. 76ff] 


Kl er. HKg,d).dt (> 0; [p] = 1jsec; [d@, P)] = sec- grd) . . . (W), 


bezeichnet werden. Wird darauf achtgegeben, daß voraussetzungsgemäß %e, 0) für alle in 


Frage kommenden Werte von o gleich Null sein soll, so lautet die aus Gl. (4) abzuleitende trans- 
formierte Differentialgleichung: 


2%v,p) 1 dw; p) — | = ae! m 
do? (#0, do +00, ;P) eg: rap /p]=sec ... (2) 
und hierin bedeutet 
„ig. pPR=y-» mit Sp] = 1 under u ee (2a). 


Diese nunmehr gewöhnliche Differentialgleichung zweiter Ordnung ist gemäß den beiden Grenz- 
bedingungen (la, b) so zu integrieren, daß 


%v,;p) = d,|p VERSATEL (26), 
%v;; p) —— d,/p Ce re a EN ern ER re Fo Beer (2B) 
ist. 
Die vollständige Lösung der inhomogenen Differentialgleichung (2) hat nun die Form 
Mesp)= 212/014 B2J2(0) 0,0. Ad (3). 


Hierin sind A und B die beiden vorläufig noch unbekannten Integrationskonstanten, während 
s_1,,(d,) die unter dem Namen der LommEischen Funktion bekannte partikulare Lösung der 
inhomogenen Gleichung (2) ist ([5], S. 234ff.). Sie kann durch die absolut und unbeschränkt 
konvergente Potenzreihe 


1 z 2 
su), 1! Due er r (22 — »2) (4? — »2) e* | 


n/2 


e RT ala 
7». sin (m v/2) Ara ya)enlaeden 


(EEE), 


definiert werden oder sonst durch die Integraldarstellung [5, S. 310 (11)]: 


TE 


. j® (v9) cos &-sno)rda na (4b). 


0 


1 
N  y.$inTv 


Statt der Funktion s_ı,,(z) läßt sich auch [(5), S. 176] die durch die Gleichung 


n]2 
y-sinTv 


S-1,(2) = s-1,(@) + 2) + I )] 


f Sin(wt 
er lz) Sn ; Iets Bora EÜBEE v0) re | e-2-Sint, Sin@d dt 
0 


y 


definierte Funktion $_1,,(2) benutzen, die ebenfalls als LommEısche F unktion bezeichnet wird. 
Im Gegensatz zu der Funktion s_ı,,(2) dieser Art ist sie so definiert, daß sie ein besonders ein- 
faches Verhalten für z2— oo zeigt. In der Tat ist dort 


By (2 _—y2) (42 — 92 
1 eh | in 8 


jim 6) 


arc (z)| < :) RE. 


16* 


RR ED Zn 
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Strebt v— 0, so gilt im Hinblick auf Gl. (5) die Darstellung: | 


x 


= [era naein Ta 
=o 
. (5b). 


2 1 
x [me -Fa+P 47-5 Fan] (megi<a) 


Auch in diesem Falle läßt sich jetzt gemäß Gl. (4a) die Funktion s—ı,, nicht verwenden. Wegen 
der Entstehung der Gl. (5b) wird auf den Text zu den Gl. 3 (12). . 3 (17) und insbesondere auf 
Gl. 3 (16) verwiesen. Der Zugang zu dieser Gleichung erfolgt dort in anderer Weise als bei G. N. 
Warson [5]. 10) 

Die in Gl. (3) noch unbekannten Konstanten A und B lassen sich mit Hilfe der beiden 
Grenzbedingungen (2 a, 8) stets eindeutig berechnen, denn es wird im allgemeinen die zu dem 
Gleichungssystem 


A J,w) + B- J_,W.) = dalp + QalPp - s-ı1,»(o) 
A. J,(v,) + B- J_v) = dulp + Oslp - s-ı,»(v) 
für A und B gehörige zweireihige Koeffizientendeterminante 


10) Is | = (50) Vi) 


A,(a, b; Yply?) = 


 sinmr (Ilm) Il) Le ren 
ee NE ra. | re.) b-vop2- 1 + E25 Cote): an 0) 
x |\L,f) Kl) s : i b2 b 


(b>a) 


nur an isoliert liegenden, endlich oder unendlich vielen Stellen verschwinden. Die Gl. (7) zeigt 
im besonderen, daß die Determinante 


a) no] 
J,(v) Yo) | Ji.(0) Y; „) 


In (b/a) 


Il 


n 


—— FE vu, =ib.yp/y? 
al wie TEN de ü BR); 
2m | cosugp- Iulbipl V +25 Coip)-dg Er ee 
0 
sowohl für positiv als auch für negativ reelle Werte von p stets selbst reellwertig ist, da I, (+ ix) 
— hl) Mir’ Imz) —0. 
Auf diese Weise erhält man schließlich unter Preisgabe der Kürzungen die folgende voll- 
ständige Formel: 
(io-Ypl®)- 


JulYpp®-ia)  Yı(Yph-ia)) 

JuYph*- ib) Yı.(Yply?-ib) 
Juli a» Yply?) pt KA Yply?) du + 95: S-ı, in (i U Yply?) 

= Julie Ypl?) Yı,lio-Ypy?) 0+9z-S_1:.(io- Ypl%) iS 
lid yo) Ynlid- ph) + By: S-rnlid- VppA)| 


in 


Die Gl. (5) läßt sofort erkennen, daß in der dritten Spalte der dreireihigen Determinante in dieser 
Gleichung an die Stelle der Funktion S__, „auch die Funktion s. -1,, von Gl. (4a, b) treten könnte, 
ohne daß dadurch eine Wertänderung der rechten Gleichungsseite einträte. 

Ist wegen I=0 nach Gl. 1(3b) auch @ = 0, so bekommt der Parameter aller Zylinder- 
funktionen in Gl. (9) den Wert 0, und überdies reduziert sich die dritte Spalte in der dreireihigen 
Determinante von Gl. (9) auf die drei Elemente 9,, 0, d,. Die Temperaturverteilung wird jetzt 
allein durch die beiden Randtemperaturen beherrscht. Umgekehrt bestimmt nach dem gleich- 


zeitigen Verschwinden von d,, d, nur noch die Wärmeentwicklung durch den Strom I den Ver- 
lauf und die Größe des Temperaturfeldes. 


N 


ER EN U TENEARTIETENUMEEN mac 5 nn 


3 


EEE 2 EIER ar A BEE a . uud u u BD 4 ml a ai le un ee» 
= A E i 
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3. Die Temperaturverteilung nach Ort und Zeit 
Um nun die eigentlich gesuchte Temperaturverteilung nach Ort und Zeit zu erhalten, ist 
auf die Gl. 2(9) die Umkehrformel der LArLAcE-Transformation anzuwenden. Ihr zufolge ist 
dann zunächst [8, S. 144ff]: 
+ioo +y 


%6) =5— er dl0o. Yplye)+ap Nm > OO... (1). 
RS 

. Auf Grund des Aufbaues von Gl. 2(9) liegt auf jeden Fall ein einfacher Pol des Integranden 
von Gl. (1) an der Stelle p = 0 vor. Er liefert bekanntlich die stationäre Endtemperatur, die 
von der Zeit nicht mehr abhängt. Um diese kennenzulernen, würde es im Hinblick auf die Gl. (9) 
genügen, die Werte der beiden Determinanten für p = 0 zu kennen, und da sich die dreireihige 
Determinante darin jederzeit durch drei zweireihige darstellen läßt, braucht lediglich der Wert 
einer zweireihigen nach Art der auf der linken Gleichungsseite stehenden Determinante bekannt 
zu sein. Hier lehrt aber ein Blick auf die Gl. 2 (7), daß offenbar die einfache Grenzwertgleichung 


3[&- 0 ypp») v(5- yo) 32m. Sin @- In bla) 
y 


I,lioypp?)  Y,lidyp®) 


+00) fürp>0 © 


besteht. Da der Grenzübergang p 0 das Argument der Funktion in Gl. (9) selbst gegen den 


Wert O0 führt, muß vorher die Funktion S_ı,, durch die Funktion s_,,, von Gl. 2(4b) ersetzt 
werden. Damit errechnet sich ohne Schwierigkeit der folgende Ausdruck für die stationäre 
Endtemperatur: 
‚sin (zu - In (o/a)) 

sin (u - In b/a) 


%o; 0) ud ‚sin (m «In (b/e)) = 


°sin(u-Imndja) ' ° 


2 sin (a In (b/o)) + sin (u-In e/a))) - 
— Oplu je Nee (3). 


An Stelle der letzten Zeile kann es auch heißen: 


Sl 


ar [0] 2% = 0) 
| sin cos (u- In Yb/a) 
Wie es sein muß, ist nach Gl. (3) 
MRS ID, ER DEE. Be ae RE E)N 
Med Le DER BE (3b). 


Das Maximum der Temperatur wird in dem Term mit dem Faktor ©, für o = Ja b erreicht. 
Es nimmt den Wert an: 


+ pl - (ee Br In Bla) — ) A RE (3%), 


95: (In(/b/a))® TÜR DENE. BOB RE (305). 
Der von der Zeit abhängige Teil des Erwärmungsvorganges wird im Hinblick auf die Gl. (1) 

und 2(9) von den Nullstellen der Determinante 
Ji,(ia-Ypy®) Yıu(lia-Yph) 
Julib- Yph®) Yiu(id- Ypl%) 
hinsichtlich p bestimmt. Solche Nullstellen sind auch für « > O0 vorhanden, denn füreinp=a- et”! 


ist dann zufolge der Gl. 2 (7) 
In (b/a) | BE S) 


Arula, b; p) = (a 


Aula,b; +iyap®)= 2m: S cosup-J(b-Yaly? -Y1 + a? — 2- a/b- Co p)- dp 


d Ima=0, 


OÖ). 
Da hierin das Argument der Funktion J, für Werte von 9=0 und 9 = In (b/a) zwischen den 
Grenzen (b—.d)- Val? und 0 schwankt und der Faktor von J, des Integranden auf dem Wege von 
der einen zur anderen Grenze bei genügend großen Werten von a verschiedenes Vorzeichen hat, 
so erhellt daraus zunächst einmal jedenfalls die Möglichkeit eines Verschwindens. 
Um auch zu zahlenmäßig verwertbaren Zusammenhängen zu kommen, läßt sich die Tat- 
sache ausnutzen, daß nach Gl. 1 (3b) u = (a - O,)"* in der Regel nur sehr klein ist. Dieser Um- 
stand legt den Gedanken nahe, etwa in Gl. (5) nach Potenzen von u zu entwickeln und beirn 
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zweiten Entwicklungsglied stehenzubleiben. Dann ist es aber unzweckmäßig, die hierbei erforder- 
lich werdenden Differentiationen am Integranden der Gl. (5) selbst vorzunehmen. Wesentlich 
vorteilhafter ist die Benutzung der folgenden Reihenentwicklung [3] für dieses Integral: Hier- 
nach ist nämlich für jedes beliebige 

Alab;tyap?) 

Ila-Yaly?) Yı,(a:Yaly?) 
Ind: Val) Yıu(b-Yaly?) 


= lyr ; >> (1,2) - 7) 


i=0 


Jin (b £ Yaly?) Ju (a- Yaly?) 
Y;.(b- Yaly?) Y;„(a-Yal%) 


J,,ı@—@)-Yah?) 
A+ 1 ER (6), 
b—a 14 n20 
(re Veh?) 
und hierin ist (i u, A) das sogen. HAnkEısche Symbol mit der Definitionsgleichung 
1 1 rg FERERIE ) 
Bee we o rG E in): De +4—ip) I, 


| r(g—A+in) 


Die Gl. (6) läßt sich in dem besonderen Falle, daß y=iu = 1/2, 3/2, 5/2... oder allgemein 
y= 1/2 + p ist, leicht kontrollieren. Die unendliche Reihe besteht dann nur aus p + 1 Gliedern. 
Der Ausdruck (i a,4) stellt also auch für u # 0 und Im (x) = 0 eine reelle Zahl dar, und für 
)— 00 zeigt er ein asymptotisches Verhalten gemäß der Relation: 


w| 


(iu) =): (2/n)!?.Cofamu-exp(A-In(Ale))-A-M2 ...... (6b). | 
Da aber andererseits für ein solches A und x > 0 
2) 
2% a: Na 
Ban) exp j (2) ee (6c) 


strebt, so ist die in Gl. (6) auftretende Reihe, weil (b — a)/2 ya b < 1 ist, absolut und gleichmäßig 
konvergent. 

Nun leitet man aber aus Gl. (6a) leicht ab, daß für sehr kleine Werte von iu=», wie 
sie nach Gl. 1(3b) tatsächlich auftreten, gesetzt werden darf: 


ra) 
=. —— ernennen: (7a), 
ee: Br wo 2) +06) ER (7b). 


Aus Gl. (6a) berechnet sich aber geradenwegs für die zweite Ableitung von (», A) nach » wegen 
Y’(1/2) = n?/2 der Ausdruck: _ i 2 S 


und so ist in erster Näherung: 


Ju(lb-Yaly) JIın(a-Yaly) 


Jo(b-Yaly?) Jo (a-Yaly*) 


Yın(b-Yaly) Yınla-yap®)| | Yold-Yap*) Yola-Yapp) | 
J ((b — a) Va/y?) 
a re A—1 1 Di ar Aa 
-H 7 e ww nd 27 4 (0, )) n r en nD “ an > 21.3: 4 
2 P2 2 ß „1 } 2yab ER NA + O(ut) 
| x (3°) Yo) 


el (0,2): N MY a) ih EA 24 
BALL (zer) 
z 5 
(8). 


werden können: Art ER 


Eu Kt, +5) u 1,04 8)) 


in) | | Y,n,rtı +9) | [5 Beier D) | 
a ern 
FE 428 (Aror -3)a +0)" j 


2yab 
Für das Symbol [3 )) in der letzten Zeile dieser Gleichung kann wiederum der Ausdruck (7b),‘ 
- (Te) gesetzt werden. Andererseits ist aber nach einer Bann Formel bei ausreichender Klein- 
heit von 6: 


= R 5< ) ©). 


en 
2 


.d 


3 is 7, ,ı@+e9) # ar ns, e 1,10) 
2 U RLEL I dy 3 
= 1 Az dena PT Ne 
Ente era) 3) Ei 
I. 2m i T,, :@ | | 
= el eg EI RD, VON ER RREL ON 
5 ) 2 (3 e) 7 ı 
2 


"und somit läßt sich der zweiten Zeile von Gl. 0) auch die Form geben: 


U) 


En 
I 3 (fo, 1 — alb)) 
ua Re u ———0. 
ö al | 2 & aan)’ 


Da nun hierin die unendliche Reihe, die neben anderen Gliedern allein den Faktor (0, A) enthält, 
die Determinante von Gl. (9) für » = 0 darstellt, so ist sie selbst der Null gleich. Für den un- 
bekannten Wert ö,, der natürlich mit der Nullstelle variiert und also von r abhängt, ergibt 
sich damit die Beziehung: 

1(fo,r (1 er a/b)) 


arte 2—1 1 bed . i+7z 

3 oa (Iarım) Sven) Sen Pe 
= nn (5 fur ap) 

er z (Ro, r (il as a/b)) 


u a ln ni a he 


ER OT WR 


(Were), 
ne 


to] co 


.® 5 2u43 
Bo e al en) RD) 
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Schließlich brauchen wir noch Formeln, die analog zur Gl. (8) die Berechnung der Funk- 
tion S_,, auf die der Funktion S_-ı,0 und eines zusätzlichen Korrekturgliedes zurückführen. Soll 
aber dabei wie bisher der Näherungsausdruck für S_ı,,(z) ein Glied der Ordnung »* enthalten, 
so ist es, wie die schon früher angegebene Gl. 2 (5) erkennen läßt, nicht zu umgehen, die Näherun- 
gen für die beiden in dieser Gleichung auftretenden Ausdrücke 


(B) s_ı,(@) und n/(2v- sinzv)- [J,l2) + I i21, 


zunächst bis zur Ordnung »* vorzutreiben. 
Nun ergibt sich nach dem Lehrsatz von TAYLoR auf Grund einer zwar recht langwierigen, 
aber sonst durchaus elementaren Rechnung für ein a < 1 die folgende Näherungsformel: 


ut Hit pan 1—a- P+ So. (RP) + a. (Pr +3 Py— P) 


1 
Hg at BP HA + pe per — Pr) + 0069) . . (12a). 
(«<1,Y=W(1 +4) ebenso für %, 9”, %” u.a.mit}=0,1,2...) 
Es wird dann bei derselben Bedeutung von Y für a =v/2 
v v y2 
(ef ey: ri er 2))= dyan: ze y 


y4 


a N De 2 13a 


und füra =» 


© (1. (z2/4)% 2 
Id) + I nn 2.|ı +5 (7 — In 22? — 9) 


an : ((#° — In z/2)? — 6 9. (P — In zj2)? +4 9” . (# — In z]2)) + 068) (13b). 


Da ferner bis zur einschließlich vierten Ordnung in » bei genügend kleinen Werten von » 


rı/2 1 mp8 7 
OR) I 
n/2 es may: 7 R : 
et + 069) a E (138) 


ist, so ergeben sich schließlich im ganzen für die obigen beiden, unter (B ä j 
! erwäh 
der Funktion S_,,, die folgenden Näherungsausdrücke: 1 REED 


En a Id) + Ja] 7 e = ze ne =. 1" + e— In 2/2)? — X’ + 72/3) 
2 . (GP — In 2/2)? — 6 9. (# — In 2/2)? +4 9”. (# — In z/2)) 
4 2 + 2m: (PP PN) +.. | ‚2 007 Wie (14) 
und 
u An Fr p3, en 


i=0 I |1 Be v\ 
: (142 Tl tar 


1. Se, wm 
DEN? 14) Na BR 
„ 2 AU | (6 Tal 


1 3 1 
RT I ; 7 ni 
4 aa Te +50) € 00) 


H. BuchHorz, Erwärmung in metallischen Kreislochscheiben bei radialem Heizstrom ... 237 


Werden nunmehr die Gl. (15) und (16) gemäß der Gl. 2 (5) zusammengesetzt, so erhält man für 
' die Funktion S_1,,(z), bei sehr kleinen Werten des Parameters, den folgenden Näherungsausdruck: 


oo 


2 
ee 06%) 


i=0 


MeTENS ee 1 3 ea 
= DI Ta lt a ee a EEE FIFT 

x A pw rm. w —_(# — In z)2)* 
+28 —n) CP — In 22? —A Pr. (Pin zj2)] +00) NE, (16). 


Auch in diesen Gleichungen ist durchweg der Kürze halber hinter den Funktionszeichen X, W”, 
Y2 usw. stets das Argument in der Größe 144 mit}=0,1,2... fortgelassen worden. Für 
die zahlenmäßige Berechnung sei daran erinnert, daß 


er aha rtn RR) (1) 


ist, und hierin bedeutet ö(n + 1,2) die verallgemeinerte Zeta-Funktion von RIEMANN. Die 


Funktion Y und ihre Ableitungen findet man u.a. vertafelt in dem Tafelwerk [7]. Aus Gl. (16) 
lassen sich zugleich die Auflösungsformeln für die beiden Integrale ? 


[0,0] [e,°} 
f pr 2. Cie. f. dt und gi e-: Eint.B.dt 
0 0 


ablesen, die anscheinend bisher nicht ausgewertet worden sind. 


Um die Beziehung für den zweiten veränderlichen Temperaturanteil zu bekommen, ist 
es danach nur noch erforderlich, die auf der linken Seite von Gl. 2 (9) stehende zweireihige Deter- 
minante mit der Determinante 3(9) in Zusammenhang zu bringen. 


Wird vorübergehend in der p-Ebene die der Nullstelle A,,, entsprechende Nullstelle mit p, 
bezeichnet, so ist demnach 


Bere, el ee (18a), 
DER ELLE ET (18h), 
und in der Umgebung von p = p, + e ist also 
ib. Y(p, + )/y? = hr +e- 1/2: h,,/p,. 


Eine leichte Rechnung zeigt dann, daß für die Determinante selbst in der Umgebung von p- 


ug) ls], Alta 


J,(ib-Yp2)  Y,(ib-Yph®) Tan en 


In (b/a) 


Jı(h,,:YI + @]® — 2: a/b.Coj y) a2 BR 
srl Zucker‘ 2 ke Lues X 19), 
27 Pr ik: 1/2: h,,-Y1 + @]b?—.2ajb- ar; |! ue- Br 


0 


geschrieben werden darf, und im Hinblick auf die letzte Gleichungszeile ist hierin das mittlere 
Gleichungsglied trotz der rein imaginären Werte von v eine reelle Größe. Für den zeitlich ver- 
änderlichen Temperaturanteil #;,,ns(0, £) berechnet sich dann damit der folgende endgültige 
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: Ausdruck: 
© J, (h, r) J, je 
drang = T' Be 
k dahe 7 5 h,, ) vs Jh, n) 


h,, ) 
er tl 9) 


e 1 
# k 4 : = Rn N) Y | b I 
x (1 + BD & v, )) Te) Y,(h,, a; De 


a a | 
ef 5 N) Br ir hr) | 


16) Hl 


Alle hierin auftretenden Determinanten können etwa nach Gl. 3(9) numerisch berechnet werden. 


Der Einfluß des von Null verschiedenen Wertes von » und des von h,,, um ö, verschiedenen 
Wertes von A,,, auf die Zahlenwerte der Temperatur kann nach den Gl. 3 (8, 11 u. 16) berechnet 
werden. 


+(% + 05: S-1,,(h,,,)) . (20). 
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Das Austauschverfahren 
der linearen Tschebyscheff-Approximation 
bei nicht erfüllter Haarscher Bedingung 


Von L. BITTNER 


Im vorliegenden Artikel wird das Problem der linearen Tschebysche [f-Approximation als lineares 
Programmierungsproblem behandelt, um auch den Ausnahmefall, daß die Haarsche Bedingung nicht 
erfüllt ist, zu erfassen. Außer einer Herleitung des „Austauschverfahrens‘‘ wendet sich diese Arbeit der 
Aufgabe zu, die Gesamtheit aller Lösungen des T sche byscheff-Problems zu errechnen. 


The problem of linear Ohebyshev approximation is here treated as a linear programming problem, 
the intention being to include the case where the Haar condition is not satisfied. Besides giving a deduction 


of In “exchange method” the paper considers the problem of caleulating all solutions of the Ohebyshev 
problem. 


B macronmeii crarbe paccMmarpuBaerca IpoÖsema mHeinol aunporcHManumm No Ye6pıI- 
IIeBy Kak IIPoÖIeMa JIHHeÜHOrO IPOTPAMMHPOBAHHA, YTOÖBL OXBATUTL TOT CHEeIHAJIbBHBbLÜ 
CtyUall, KOTIA YCIoBHe Naapa He BbINOAHeHO. KPoMe BhIBoA1a „‚MeTona 3aMmeHpI“ aBTop 32- 
NAOTCH IICIIBIO YCTAHOBUTL COBOKYIIHOCTB BCexX PelmeHnH ITPoOIeMbI YeösimeBa. 
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1. Einleitung 

In [1] publiziert E. STIEFEL das von ihm entwickelte „Austauschverfahren“ zur schritt- 
weisen Lösung von Aufgaben der linearen TSCHEBYSCHEFF-Approximation. Herleitung und 
Wortlaut der Rechenregeln sind jedoch zunächst an die Voraussetzung gebunden, daß die homo- 
genen Teile von je m linearen Funktionen (m = Anzahl der Unbekannten) linear unabhängig 
sind (HaArsche Bedingung). In einer darauffolgenden Arbeit [2] stellt STIEFEL fest, daß sein 
Austauschverfahren in einem gewissen Sinne der dualen Simplexmethode, angewendet auf ein 
zugeordnetes lineares Programmierungsproblem, äquivalent ist. Diese „Äquivalenz‘ soll hier 
ausführlicher als in [2] betrachtet werden, um auch alle die Fälle zu erfassen, in denen die oben 
erwähnte Voraussetzung nicht erfüllt ist, Bei nicht erfüllter HaArscher Bedingung können analog 
zu ‚den entarteten zulässigen Basislösungen eines linearen Programmierungsproblems auch 
„charakteristische Relationen“ mit einigen Koeffizienten A, = 0 und entsprechend den Basis- 
zyklen Wiederholungen schon benutzter Referenzen vorkommen, wenn man nicht besondere 
Regeln zur Bestimmung der aus der Referenz auszuscheidenden Zeile beobachtet. 

Die Formeln (42.1)—(42.7) und (51) enthalten das Hauptergebnis der vorliegenden Arbeit, 
nämlich die modifizierten „Austauschregeln“. Abschnitt 5 weist auf die tabellenmäßige Fassung 
des Verfahrens, die ‚reduzierten Tableaus‘‘ von STIEFEL, hin. Im Abschnitt 7 werden Wege 
angegeben, wie man zu einer Darstellung der Lösungsgesamtheit der TSCHEBYSCHEFF-Aufgabe 
gelangen kann. 


2. Problemstellung 


2.1. Es handelt sich also um das Problem, zu dem unverträglichen, reellen Gleichungs- 
system 


NZ; +G—=O ee I I HE 
k=1 
ein System von Variablen x,, 2, ... , 2, zu finden, für welches 
ln, a nn = Max m) ne el 
ein Minimum wird. Wie in [2] und [3] dargelegt wird, ist diese Aufgabe dem folgenden Problem 
äquivalent: 
2.2. Es ist das Minimum der linearen Funktion 


I ar ae ee: are a) 
unter Berücksichtigung der Nebenbedingungen 


m m 
Ni nz Im+1 = 2& Ad;r X ns Lm+1 4 C; = 0 D —Mi 4 Ra & Air X + Im+1 € Pi 0 (4) 
k=1 ki 


((=1,2,3,...,m) aufzusuchen. Diese Aufgabe stellt jedoch noch kein typisches lineares Pro- 
grammierungsproblem dar, weil die Variablen x,,%,...,2%„ auch negative Zahlenwerte an- 
nehmen dürfen. Die Überführung in ein solches Programmierungsproblem kann natürlich auf 
sehr viele Weisen durch Einführung neuer Variablen vorgenommen werden. Im folgenden wird 
eine ganz spezielle Möglichkeit gewählt. Zunächst notieren wir einige Abkürzungen. 

2.3. Es sei W = (a;,), r = RM) der Rang von Y\, WU, eine r-zeilige Teilmatrix von W des 
Bangesırı, 2. die (n 7 r)-zeilige, Restmatrix von Wir = (ZT, a... , La) C = (Cs Es Jans a)» 
der den r Zeilen der Matrix U, zugeordnete Teilvektor von c, c, der denn —r Zeilen von W, ent- 
sprechende Teilvektor von c. Offenbar giltr <m,r <n. Jede Zeile a” von W, kann als Linear- 
kombination der Zeilen von W, dargestellt werden: a® = t’-W,; ’ bedeutet hierbei den Zeilen- 
vektor der r Linearkombinationskoeffizienten. Infolgedessen gilt 


7, Fellbaına.. Breit 
wenn man unter X die aus den (n — r) Vektoren t’ zusammengestellte Matrix versteht. 
Es bedeute ferner e; = (1,1,...,1)’ (i Komponenten), 0; = (0,0,..., 0) (i Komponenten) 
a,. den i-ten Zeilenvektor der Matrix W, G; eine i-reihige Einheitsmatrix. 
Durch 
4 yı =S U, t u er . Cm+1 -F (ı . . . . . . . . . . . . . (6) 


wird ein r-zeiliger Vektor definiert. Dieser Vektor ist nicht-negativ, wenn die Variablen x,, 
X + +, %n+ı dem ersten Teil der Nebenbedingungen (4) genügen; er nimmt jedoch jeden mög- 
lichen Vektor des r-dimensionalen Raumes an, wenn die Variablen x,, X, . . . &n-+1 frei veränder- 
liche Unbestimmte sind. Es gilt wegen (5) und (6) 
U; c in en—r ' Im+1 IF (5 = T Yı -r- (ee Zi T &) Er Im+1 -r (a = T Cı » | 
— A, t + er, T Im-+ti rue Cı — er: Yı En 2 er 7 Im+1 P) . . (7). 
— WU, % = Een—r * Im+1 (2 == 2 T Yı -F (643% + T e,) * Im+1 r T Cı <>; (a | 


* 
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i ten von y, und der Variablen x.+1 
. Aus (4) folgt noch, daß x > 0 sein muß. In den Komponen 19ı 
hat nn ae Erfüllisein yon (4) gerade r+ 1 nicht-negative Variable, mit deren Hilfe sich 
das Problem 2.2. in ein typisches lineares Programmierungsproblem umformen läßt. | 
2.4. Äquivalente Fassung der ursprünglichen Aufgabe: Es ist das Minimum der linearen 
Re egal. an al. Aa Rn 


unter Berücksichtigung der Nebenbedingungen 


T (ee a: T er yı T Cı == Ca i 
u, 2.0, [| = 0, & 7? RE ö (9) 
| T En—r -H T er Im-+1 (Ca ——z T Cı 


aufzusuchen. Die Äquivalenz ist im folgenden Sinn zu verstehen: Jeder Lösung der Aufgabe 
2.1. bzw. 2.2. entspricht gemäß der Substitution (6) genau eine Lösung der Aufgabe 2.4.; um- 
gekehrt aber wird jeder Lösung von 2.4. im Falle r < m durch (6) eine 00” "-Mannigfaltigkeit 
von Lösungen der Aufgabe 2.2. bzw. 2.1. zugeordnet. 


3. Grundtatsachen des linearen Programmierens 


3.1. B sei eine endliche Matrix, p und ) seien Vektoren mit soviel Komponenten, wie B 
Spalten besitzt, x und 3 dagegen Vektoren, deren Komponentenanzahl mit der Zahl der Zeilen 
von ® übereinstimmt. Bekanntlich haben das primale Problem 


p’-y = Min! für, 39T y>D, Ve 
und das duale Problem 
v3 Marl Sun Seh, s>D, (3 = 7), wer 


— wenn überhaupt, dann — gleiche Extrema. Durch Definition von nicht-negativen Schlupf- 
variablen, die wir zum Vektor 3: = p’ — 3-3 zusammenfassen, und die Vereinbarung, daß 
3 t 28) EEE f 
ve ( ) q = (hier o = Nullvektor der Dimension von 3,), % = (<) (& = Einheitsmatrix der 

0 ) 
Dimension von 3,) bedeute, kann das Problem (11) in das erweiterte duale Programmierungs- 
problem 
q-3= Mal für 3 8=P», =D (=? De un . (12) 


übergeführt werden. Die Zeilen f; von % sind den Komponenten von 3 und q in eindeutiger Weise 
zugeordnet. Mit 5* werde eine zulässige Basis, mit 3*’” = p'($*) > v’ die zugehörige zulässige 
Basislösung und mit q* der zugeordnete Teil des Vektors q bezeichnet (siehe [4]). Die Simplex- 
methode erlaubt nun, mit Hilfe einer als vorhanden betrachteten zulässigen Basis %7 und zu- 
lässigen Basislösung 37° = p’- (31)! eine im allgemeinen ‚bessere‘ zulässige Basis 53 mit der 
zugehörigen Basislösung 32° = p’- (53)! zu konstruieren. 7, q$ mögen die den Vektoren 3} 
und 3% zugeordneten Teilvektoren von q sein. 
3.2. Ein Schritt der Simplexmethode besteht aus der 


Feststellung a), welche Zeile f, von % als Zeile in die neue Basis aufgenommen werden soll. 
Man bildet zu diesem Zweck den Vektor 


TEN) ER 
und stellt fest, welche Komponente fu F)'q7 —q< 0 ist; die zugehörige Zeile fr bzw. — 
beim Auftreten mehrerer negativer Komponenten — eine der zugehörigen Zeilen ist als Zeile 


der neuen Basis zu verwenden. Wenn alle Komponenten des Vektors (13): > 0 sind, stellt 37 
Bes Aine; optimale Basislösung dar. Die Rechnung wird bis zum Eintreten der Relation 
Re & ) 0” — q = 0 fortgeführt, obgleich eventuell schon vorher eine entartete Basislösung 3* 
optimal sein kann, ohne daß alle Komponenten von (13): > 0 sind). 

I, 


Feststellung b), welche Zeile f,, aus der Basis $+ — fi, — | zu entfernen ist. Man bildet 


zu dem Zweck die zulässige Basislösung zf’ = pP’ mit den Komponenten z* und den Zeilen- 
!) Das ist der neg 
gehörigen Variablen. 
?) Daß die Bedingung (3 


al 7 * > : N) * * 1x . . 
ta et, 1 q q= 0 im Entartungsfall nur ein hinlängliches Kriteri ü i 
malität darstellt, ist in [4] anscheinend verkannt worden R z BU SUr IET 


ative Vektor der Koeffizienten der Zielfunktion, ausgedrückt durch die nicht zur Basis 


» 
i 
h 
| 
ö 
| 


asl TEE SZ EZ Ze j W 
= u ’ 
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vektor d; = fr: (31)! mit den Komponenten D;,, und stellt fest, für welchen Index i = I! das , 


min a TREE 
(1:%,5>0) Vi,r 


angenommen wird. (Falls alle v;,, = 0 sind, existiert für die Aufgabe (12) kein eigentliches 


| Maximum; die obere Grenze ist oo. Dieser Fall kann in der-sogleich zu betrachtenden Aufgabe 


wegen des Vorhandenseins eines eigentlichen Extremums nicht eintreten). Die zu dem Index I 
gehörige Zeile f,, bzw. — wenn das min mehrfach angenommen wird — eine zugehörige Zeile ii 
ist aus 1 zu streichen; die Restmatrix ergibt zusammen mit der in a) ermittelten Zeile }; die 


neue Basis %. 
Beim Übergang von 37 zu 3% nimmt die zu maximierende lineare Funktion (12) um 


2. 2-1 4-19) +0 et — dr g) 


x >0. für 9=#0 
ge rar RAN N, 
(1: — Si) gr) < 0 Für ee (15) 


3.2.1. Um die Wiederkehr schon benutzter Basismatrizen zu vermeiden, darf man die aus- 
zuscheidende Zeile bei mehrfacher Annahme des min (14), wenn dieses min = 0 ist, nicht will- 
kürlich unter den anderen möglichen Zeilen auswählen. Bei mehrfacher Annahme des min 0 
wollen wir von einer „ernsthaften Entartung‘‘ sprechen. Man geht beim Auftreten einer erst- 
maligen „ernsthaften Entartung‘‘ zu einem abgeänderten Problem, einem e-Problem über, 
das dadurch ausgezeichnet ist, nur Basismatrizen zu haben, die auch Basismatrizen des ursprüng- 


ZU. 


‚lichen Problems (12) sind, und zugehörige zulässige Basislösungen zu haben, die nicht-entartet 


sind. Sei %0 eine zulässige Basis von (12), für die eine lokal erstmalige ernsthafte Entartung 
eintritt; unmittelbar vorher möge das min (14): = 0 gefunden oder — falls es = 0 ist — nur 
einmal angenommen worden sein. %, sei die umgeordnete Matrix %, in der die ersten Zeilen 
in irgendeiner Reihenfolge die Zeilen von %5 sind. Es sei & eine sehr kleine positive unbestimmte 
Zahl und € = (8, &, &®...) ein Zeilenvektor aus den aufeinanderfolgenden e-Potenzen mit soviel 
Komponenten, wie % Zeilen hat. Dann lautet das e-Problem 

g’.3()=Maxl für zZ) S=P+E- (ed) SD GE?) - - (1). 
Auf das e-Problem (16) ist die Simplexmethode genauso anwendbar. Zur Ermittlung der aus der 
Basis %f zu entfernenden Zeile f,, hat man die zulässige Basislösung zt’(e) = p’- (Fr)! + EFF)" 
mit den Komponenten z?(e) zu bilden und festzustellen, für welchen Index i = I! das 


* 
ED DL er ee Sea. I 
(1:%,5>0) ‚Pi,r 

angenommen wird. Es gibt jetzt nur einen einzigen Index I. 6(e) ist > 0. Die erhaltene Basis 
5% ist auch eine Basis von (12), und inbezug auf das e-Problem hat die zu maximierende Funk- 
tion um 6(e) (9. — fr(&7) '- qi) zugenommen. Wir ermitteln die weiteren Basismatrizen so- 
lange mit Hilfe des e-Problems, bis wieder einmal ein von 0 verschiedenes min (14) bzw. ein 
min (17) mit 0(0) = 0 auftritt oder die Rechnung abgebrochen werden kann, weil der Vektor 
(13): vo ist. Danach gibt man das e-Problem auf, weil die zusätzlichen Rechnungen zur Be- 
stimmung des min (17) unbequem sind, und definiert später notfalls wieder ein neues e-Problem, 
wenn nochmals eine lokal erstmalige ernsthafte Entartung auftreten sollte. Weil entweder 
die lineare Funktion (12), die Zielfunktion, monoton zunimmt oder die lineare Funktion (16) 
wenigstens um Größen der Ordnung e* (h 1) zunimmt, während die lineare Funktion (12) 

konstant bleibt, können schon benutzte Basismatrizen nicht mehr wiederkehren. 
3.3. Mit Hilfe einer optimalen Basis %* des dualen Problems, welche die Bedingung 
SSHTga*—q>o erfüllt, gewinnt man bekanntlich eine optimale Lösung des primalen 


Problems durch N a er er BR an ae er etyt 


4. Das Tsehebyscheff-Problem als duale Programmierungsaufgabe 


4.1. Wir wollen mit 2,, 23 . - . » 2, 22+1 die (r + 1) Schlupfvariablen des zu 2.4. gehörigen 


erweiterten dualen Problems, mit 2,41; - - » » Zus 2n4% + + » Z2n+ı dagegen die übrigen Variablen 
bezeichnen. Dieses erweiterte duale Problem lautet nun zunächst - 
D, 
Tu — 0% 
a ae 0 N ee I) 
um 
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mit den Nebenbedingungen 


EC 0, 
T En—r T er 
1.5 =34- D, 1 TE N = Dan+ı - - (20). 
—€&, 2 } 
—tT YottTe 
6, 
Is | 
Weil eine zulässige Lösung 5 also dem Gleichungssystem 3» o,| = vo, genügen muß und 
—E 
infolgedessen | —tT 
6, [ G 
T “ c 
a Duni tarr“ (ı=-# 0 
—E, un.) uf 
J u rs 
sowie ö % 6 
En—r — T er | Een—r z 
5 1 = pt e: dr | = 3: eantı 
Dre; | & —& 
nr tTe, eu: —{T 
gilt, läßt sich die betrachtete Aufgabe auch in der symmetrischen Form schreiben 
2 
5 .| 0)=Mal für 5 Ir e&ontı]l = (0 1), 52 Dan+ı - (20). 


c Tys (67 
\- 2 
Unter den Zeichen %, p und q sind jedoch im folgenden stets die aus (19) abzulesenden Matrizen 
und Vektoren zu verstehen (vergleiche mit (12))! 


D, a 
3” -() aut 7=2, 23=2,..,2%=%, 24-4 


ist augenscheinlich eine mögliche Basislösung. Eine Basislösung hat hier immer r + 1 Kompo- 
nenten, eine Basis ist eine (r + 1,r + 1)-Matrix. Es soll noch ein für allemal vermerkt werden: 
Wenn eine Matrix bzw. ein Vektor mit einem * versehen ist, dann ist die durch die Basis ur 
bestimmte Auswahl von (r + 1) Zeilen bzw. r + 1 Komponenten gemeint. 

4.2. Um künftig Wiederholungen zu vermeiden, wollen wir zunächst eine mehrfach ge- 
brauchte Identität notieren. Es ist 


Er eL, 
‚ T Gr — Te 
Jı re PER UEt+ Er imr + Ci 
ll = 0, 1 : 
“ G Im 
Im+1 —E 26% or 
—T En—ı T er 
U, 1% + e,* Y-} 1 + qı G, 
T A, c + Enrar 8 Im-+1 -)- T Cı ü U T 
= Tm+1 = Om|E-+ Eantı"Imtıt Or | + cı (21). 
To U, ” m Er * Um-+1 — cı — A 7 02 
—T Ar -+ En—r * Im- ı— Te —T 


Entsprechend ist 
AU\* 


T 

= [sh \=| «| rer ee 
eo 
=; 


Km Fi L +1 A 


7 Ei) el a SD ze De 
- 
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R 4.3. Um zu einer Deutung einer zulässigen Lösung, insbesondere einer zulässigen Basis- 
lösung zu gelangen, multiplizieren wir das System von Nebenbedingungen (20) mit dem Vektor 


1 . s 2 ; , 
® ). in welchem wir die Komponenten %p %y..., %m+ı als frei veränderliche Unbekannte 
m-+1 3 


betrachten. Wir finden mit Rücksicht auf (21) und (22) 


A 
5 )=r DL eo In ed 
— — A 
bzw. \ - 7 g\%* wa (23). 
Borken | Im) E43" er ihm 
“m-+1 gr 


Hieraus gewinnt man durch Koeffizientenvergleich in bezug auf x,, &,, ..., &n+ı (und eine weiter 
unten angestellte Überlegung) das 


1. Ergebnis: Eine zulässige Lösung 3 ist genau eine durch 3° - ea, = 1 normierte, nicht- 
negative Lösung des homogenen Gleichungssystems 
A 
% ® De = Dir DEISPEOr BORD, re (24); 
— U 
/ 
eine zulässige Basislösung 3* ist genau eine durch 3*- e,.ı = 1 normierte, nicht-negative 
Lösung des homogenen Gleichungssystems 
NERHNT 


| an u en a rd nee he 
4 


Daß auch umgekehrt jede durch 3 egn+ı = 1 normierte, nicht-negative Lösung von (24) 
dem System von Nebenbedingungen (20) genügt, erkennt man so: In WU, gibt es eine nicht- 
singuläre (r, r)-Matrix W,ı- Die aus denselben Spalten wie W,, bestehende Teilmatrix von W, 
heiße W,,. Offenbar ist A, = TW,ı- Aus dem Bestehen von (24) darf man somit folgern 


Azı T T 
er ee res. 
>= HU =; G, un: (07 


Zu Au) Ic£ T & T 
Damit ist das System von Nebenbedingungen in (20’) erfüllt und deshalb auch in (20). Gleiches 
gilt für eine zulässige Basislösung, die ja eine zulässige Lösung mit höchstens r + 1 von O ver- 
schiedenen Komponenten darstellt. Wegen der Nicht-Singularität von 5* ist 3* eindeutig be- 
stimmt; daher muß auch die Matrix des Gleichungssystems (24°) im Verein mit der Normierungs- 
vorschrift 4*’ - e,+1 = 1 nicht-singulär sein. Daraus läßt sich eine charakteristische Eigenschaft 
der Systemmatrix (24°) erschließen. 
\* 
Zum 1. Ergebnis: Die Matrix 0„,| hat den Rangr, die Zeilen sind also linear abhängig, 
— A 
und die von Null verschiedenen Linearkombinationskoeffizienten haben gleiches Vorzeichen. Nach 
A\* 
STIEFEL möge eine solche Auswahl von (r + 1) Zeilen, wie sie in der Matrix 0 | vereinigt 
en 

ist, eine Referenz heißen. 

4.4. Wie in 3.2 dargelegt worden ist, verlangt die Simplexmethode bei Kenntnis einer 
Basis %* die Berechnung des Vektors 3(3*) !q* —g. Wir wollen nun den Vektor A ae 
yı 


) darstellen. Es sei also jetzt 
Lm-+1 


als speziellen Vektorwert von ( 


Mel ©, Din - 


Km +1 


(FHT:gr = | . (25). 
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i j j ist, daß wenigstens ein Größen- 

Die Variablen sind nunmehr durch (25) festgelegt. Klar ist, or 
system X %y +.» &m+ı existiert, welches (25) erfüllt; &m+ı wird eindeutig bestimmt, m Rn r 
Komponenten von x dagegen bleiben willkürlich. Die Variablen x,,..., %„;ı müssen der Be- 


ziehung 
a (ar 
T Cı ae Ca au i T 
= 0 = Om] :E+ Errı Iimsı rt d, | GG 
0, u iA 
(a Betz T Cı je m z ) 
gehorchen, d. h. es muß sein ; 
DJ * c * 
Dr -c+ 0 + &41 m = 0 u ee . (26). 
— A —() - 
Nach Multiplikation mit 3*’ folgt wegen (24°) 
BR c Pr 
Inte Ol... ee a 
c 


Die durch (27) definierte Größe heiße wie bei STIEFEL „Referenzabweichung‘; sie stellt gleich- 
zeitig den Wert der linearen Funktion (19) dar. Die Komponenten des Vektors r sind danach 
aus (26) zu errechnen. Führt man noch die Abkürzung 


Nele. 5 en Be 
ein, so läßt sich die Bedingung (26) in der Form 


Y) * C * 
0 ei ee 0 Ent: sun 2 a Bu ae (29) 
—) —t 
schreiben. Die zur betrachteten Basis gehörigen linearen Funktionen + n; müssen für die ge- 
suchten Größen X, . . . , %, denselben Funktionswert haben. Ein Vektor r, der (26) erfüllt, möge 


nach STIEFEL ‚„Referenzzentrum‘ heißen. Indem wir weiterhin unter r, 2&„+ı ein spezielles 
Lösungssystem von (25) verstehen, können wir folgendermaßen weiterschließen: 


EB 


AMıtemarte a z 
Er) q* >= > D) | “ DE ) = Om a O0) + ean4ı In+1 
m-+L1 er er 
9 
= 0 + Eon Dial) ir nr a en ae Kerpen ee (30). 


—) 


Die Komponenten 7, von 9) sind hierbei an einem „Referenzzentrum‘ berechnet. Man 
kann nun ohne Beschränkung der Allgemeinheit annehmen, daß der Zielfunktion (19) bzw. (20) 
durch die gerade betrachtete Basislösung 3*’ = p/(F*)! ein nicht-negativer Wert erteilt wird; 
dann ist auch 2.41 = 0 in (30). Wenn dies nicht so sein sollte, läßt es sich leicht ändern. Aus 


einer zulässigen Lösung 5 mit negativem 3’. q läßt sich mit Hilfe der Festsetzungen 2, ::= n+1+j 
(=12,3,..,7), Zn Saiı Zuat y=3(j=12,3,...,n) eine zulässige Lösung 3 mit 


positivem 3-4 = — 3’: q gewinnen. Im Verlauf der Anwendung der Simplexmethode kann 
die Zielfunktion nur wachsen oder konstant bleiben; x,,.;ı wird dann also auch nach jedem 
weiteren Simplexschritt = 0 bleiben. Wir fassen die Aussage der Formel (30) zusammen zum 

2. Ergebnis: a) Die zulässige Basislösung 3* ist sicher optimal, wenn für en. 
sowohl 7; + Zn+ı Z0 als auch — 7; + m+ı 20, d.h. || S aum+ı ist. 

P) In die neue Referenz ist eine Zeile a) oder — a; aufzunehmen, jenachdem nz + Im+ı < 0, 
also 2.41 < —nı = |n.|, oder — Ne + %mrı < 0, d.h. &u+ı < Ne = |Nel, ist. Mit anderen Worten: 
Falls |nx| > &n+ı ist, kann die Zeile — sgn nn, a, in die neue Referenz aufgenommen werden. 

Die Größen n; sind hierbei die Funktionswerte der linearen Funktionen n,(&) für ein be- 
liebiges Referenzzentrum x der alten Referenz. 

Gemäß 3.3 kann man aus einer die Bedingung a) erfüllenden optimalen Basislösung mit 
der Basis 5* eine optimale Basislösung des primalen Problems 2.4 und somit auch eine Lösung 
der ursprünglichen TSCHEBYSCHEFF-Aufgabe 2.2 herleiten. 
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Diese optimale Lösung für 2.4 lautet 


=(" )-@Ntar. 


Lm+ 
Durch Auflösung von 


B 3 r er ICm+1 1E C 


— (Zr) 1 g* 
A Erg 
nach r erhält man eine Lösung von 2.2. Dieses Gleichungssystem stellt aber gerade die mit (25) 
identische Definitionsgleichung der Referenzzentren dar. Wir kommen daher zu der folgenden 
Ergänzung des 2. Ergebnisses: Wenn a) erfüllt ist, stellen die zugehörigen Referenzzentren Lösun- 
gen der ursprünglichen TSCHEBYSCHEFF-Aufgabe und die zugehörige Referenzabweichung das 
gesuchte Minimum der maximalen Abweichung der linearen Funktionen n; dar. 

Um die aus der alten Basis bzw. alten Referenz auszuscheidende Zeile zu ermitteln, hat 
man nach 3.2. b) den Vektor 


v4 = fi . (Fr) bzw. Dnrırk = er . (Fr) A er og hs (31) 


zu bilden, je nachdem 7, < 0 bzw. n, > 0 befunden wird. 
Der gesuchte Vektor hat also der Beziehung 


=. % = 9:5” bzw. Inrı322 = nr14e 5 = Dat14n 5° Al (31’) 
zu genügen; d, repräsentiere dabei den Vektor (61,2 Ö2,1 - - - » dantı,n)'’; On+ı+. entsprechend 
A » Öan+i1,n+ı+r). Die rechtsseitige Multiplikation von (31’) mit dem Vektor 
[* Eu t .) ‚in welchem die Variablen nun wiederum als freie Veränderliche fungieren 

Im-+i1 
mögen, liefert die Identität (siehe (21), (22)) 
= 5° 
T AN= < 
ar 5 -+ Imtı + dk 9, |o=d- Om] E+ dr &r+ı &mtı + Or* | > 
>: G, — A TR F, 
u, —T 
aus der durch Koeffizientenvergleich schließlich die Gleichungen 
DIR * 
— D;- Oel; le: er +1 
— A 
d R 3 
> © Ne. A en 
T T 
d; Dee Dre €; 
ze —E, 
\_ x 
u 
folgen. Die letzte Gleichung erweist sich wegen a4=d;: D„ | als einfache Folgerung des 
— A 


ersten Teiles des Gleichungssystems (32), wie analog zu einer im Anschluß an Formel (24) ge- 
führten Überlegung gezeigt werden kann. Aus dem Bestehen von (32) kann rückwärts auch 
auf das Bestehen von (31) geschlossen werden. Die Bedingungsgleichungen von n+1+2 lauten 
a‘: 
PER A = On+1+k 5 Din ’ i= Dn-+1+% °&,+1 
— 4 
Die Komponenten z* der zulässigen Basislösung 3*, die Komponenten v;,, von dr en 
von d„+153, und die Zeilen + a; einer Referenz sind einander eindeutig zugeordnet. Wir fassen 
Ar 

zusammen 


Be 


17 
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3. Ergebnis: Aus der (alten) Referenz zu entfernen ist eine solche Zeile der Referenz, die 
der Basisvariablen zugeordnet ist, für welche das 


Rz 


im Falle 7. <0 


min 
(:9,,5>0) Pi,k 
bzw. 
z > 
min —— im Falle .>0 


(#:%;,m+1+%> 0) Pin+1+k 


angenommen wird. z#, v;, und D;n+1+x Sind aus den Gleichungen (24’), (32) und (33) zu er- 


rechnen. 

4.5. Obgleich die Rechnung nun vollständig mit den ursprünglich gegebenen Größen Ars 
c; durchgeführt werden kann, wirkt sich die Verdoppelung der Gleichungssysteme auf die Hand- 
habung des Lösungsverfahrens erschwerend aus. Diese Schwierigkeit kann man beheben. 


Um das Austauschverfahren in der Fassung von STIEFEL zu erhalten, treffen wir folgende 
Vereinbarungen: Wir gehen von einer Basis bzw. Referenz aus, die lauter verschieden indizierte 
Zeilen a,, enthält; zwei Zeilen a; und — a; dürfen also nicht vorkommen. Als Folge der An- 


wendung der Simplexmethode besteht auch jede weitere Referenz aus lauter verschieden indi- 
zierten Zeilen. Wir wollen außerdem -annehmen, daß die Zeile o,, in der Referenz nicht vor- 
kommt. Die Zeilen einer Referenz, umgeordnet nach wachsenden Indizes, ohne das eventuelle 


Vorzeichen —,, fassen wir zur Matrix 
’ 
A 
a, 
Er ae 
Yr — . IS Zu. 2, ug on . . z . . . (34) 
Er 
TTS 
zusammen. Die Komponenten einer Basislösung sowie die des Vektors dv; bzw. d„ı117; Seien 


ebenfalls so umnumeriert, daß wachsenden Komponentenindizes auch wachsende Indizes der 
zugeordneten Zeilen der Matrix W* entsprechen. Das Vorzeichen der Zeile Ay; in der Referenz 


’ gy\* 
0] heiße ß;. 
a 
Mit Hilfe einer beliebigen reellen Zahl o lassen sich Größen 
A:=0-ß-z ee 
definieren, die, zum Vektor 


% = (dis An nung Ar 1)" 0, RR Em u Kae nn (35) 
zusammengefaßt, offenbar dem homogenen Gleichungssystem 
1 UR ed 1 HU De 2 


genügen. Umgekehrt erhält man aus einer beliebigen Lösung eines Gleichungssystems vom 
Typ (36), in der A* eine beliebige aus r + 1 Zeilen bestehende Teilmatrix von W des Ranges r 
sein darf, mit Hilfe der Festsetzungen 
A| 

a | T a « om 

= oT V=l23,..5r FD 
2 A| 
je] 


Komponenten einer Basislösung, die den Zeilen I entweder mit den Vorzeichen + sgn A, oder 

En een ra sen A, für A; = 0 und beliebigem Vorzeichen für A; = 0 zuzuordnen sind. 
amit  Referenzabweie o Hm ancfä \ a. za Q nzoi r n 

ren ie elerenzabweichung positiv ausfällt, setzen wir das Vorzeichen von a,,in der Referenz 

P:=T-sen}, für #0 | 


/ ri -ı 
= son (ic) | 2 in. 2) a at ae a 
—1 


En Die a: P; für 4, —= 0 dürfen wir nach Belieben festlegen. Falls die betrachtete 
eierenz jedoch durch die vorangehende re g | X ) 
Te Re ee angehende Anwendung der Simplexmethode entstanden ist, 
liege en Zeilen zuzuordnenden Vorzeichen bereits fest. Die Vorzeichengebung (37) stimmt 
übrigens mit der von der Simplexmethode geforderten Festsetzung überein. 


mit 
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Aus den Formeln (32) und (33) ergibt sich mit Rücksicht auf unsere Vereinbarung hin- 
sichtlich der Indizierung, daß v,(0%n414.) eine durch 1% 841 = 1, (Wizisr- &+1 =1) nor- 
mierte Lösung des Gleichungssystems 


Pı =0 
er Pa ’ 
a: = dr * “ -Yr: = m, -U* 
bzw. ’ ee a RE (38) 
ßı 0 
— A = Dnr+ı+k' P. : » J* 


0 Br+i 
ist. Aus einer speziellen Lösung 1m. = (w1,% 42% - - - » Ar+ı,»)’ des Systems (38) und der all- 
gemeinen Lösung s-| (s reeller Faktor) des homogenen Gleichungssystems (36) erhält man 
demnach db; und d„+1-+z in der Form 


(m +s’-D, 
Ber Er EN 5) 


D; 


Pı 
Dn+1+ = i$ N) 


Ba 


Die speziellen Parameter-Werte s* und s- können aus der Normierungsbedingung errechnet 
werden. 

Es werde zunächst angenommen, daß die gegenwärtig betrachtete Basislösung nicht- 
entartet, also jedes A; #0(i=1,2,3,...,r +1) ist. Wir ersetzen nun zf und v; x (I,n+1+2) 
durch Größen in A, und w;,,„. Mit Rücksicht auf (37) und (39) gilt 

min min 


; et r+1 rl ; nn: 
mar re AN ABl sta) || maxß:- sen hl > a) 
j=1 j=1 in 


r+1 Br ar ; a 
IA; | max 77 +r:s 
j-1 { ee ) } “ 
v:T +st] >0 
); 
[2 
1 1 
= . I RETAIL 4 A) 
r-+1 7 k + 
1A;| max (7 En + T+S 
j=1 i=1...r+1 Ü 
und 
y* 1 1 
min =. —- ? —— .....(40). 
20; > un+1-+k is Re 
(l:d;,n+1+r>0 "und 1-+k 1A,| max Er „Ei ER: 
j-1 i=1...r+1 4; 
* n„* 
Die Ermittlung eines Index i, für den.das min — (min ee angenommen wird, ist 
i,k i,n+1+k, 


x = A B . .. Ki,k 
damit auf die Ermittlung eines Index i, für den das max Kt sen nk* ei angenommen 


V ) 
wird, zurückgeführt. Wegen der Beschränktheit der Zielfunktion (19) gibt es immer wenigstens 
ein positives v; ,. Nun werde auch zugelassen, daß einige A; = 0 sind; die zugehörigen Kompo- 


nenten v; ; (Yinyı+.) lauten —sgnng- Pi Air. Wenn sich unter diesen Größen wenigstens 
I $) a b) 2 


r 2 “ii . ii . FE . 2 
eine positive befindet, ist das min - (min Ge) —= (0) und wird für einen Index i 
D;,% Dj,n-+1-+k 
mit der Eigenschaft 4; = 0 und — sgn n + Pi Wi, > 0 angenommen, Falls die Indexmenge 
„* N „* 

7 . . “ii . 577 E , 

eig, son ne Bist >% leer ist, kann ‘dieses min — (mi —— ) nur für 
a BUN Laer ; ; D;,k Dj,n-+1+k 


nl: 
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Indizes mit A; + 0 angenommen werden und ist selbst # 0. Man erhält also ähnlich wie in (40) 
bzw. (40°): wenn I leer ist, dann ist 


* 1 1 
min Sr = +1 . 3 er. se De EN (41), 
su r>0R dirk 4 fi,k act 
(> Zu. maxi er 
Fat {1:4 +0) i 
bzw. 
4 1 1 , 
En x a Le. = — ..7.li 
fir,ntı4r> 0 Pünti4tE 5 4j| max pe 4) +r.& 
j=1 {i:7, #0) hi 


Insgesamt läßt sich ein Schritt des Austauschverfahrens aus den folgenden Rechenschritten 
aufbauen: 


1. Schritt. Bestimmung der „charakteristischen Relation“ etc: 
ı Auflösung des Gleichungssystems 


i& . + = Di (tl = (Ay, An ...3,3 Ar+1)') au ae ie (42.1), 
Berechnung der zugehörigen „Referenzabweichung‘“ 


we 
IV - c*] j= i 7 42.2 
h= m le Be Re 
zi zZ 


' Festsetzung der Vorzeichengrößen ß; für die Zeilen a,, am Anfang 
ziesgn lich), ße u sgnd) are 

Be= Hl beliebig für 4=0 ' 

oder Übernahme dieser Vorzeichengrößen aus der vorangehenden Referenz. 


und 


ei) 


2. Schritt. Ermittlung einer in die neue Referenz aufzunehmenden Zeile a;: 
Berechnung der Werte der Funktionen n; für ein gewisses Referenzzentrum, das durch 
die Gleichungen 

N, = —Birh i=1 23.5.3 PD EI 
definiert ist, 
Auffinden eines Index %, für den nr >höailt, 
Notierung der Vorzeichengröße ß := — sgn n; der Zeile a; in der neuen Referenz, 
oder Abbrechen der Rechnung, wenn In;! sh (=1,2,...,n) gefunden wird, und 
Berechnung eines Referenzzentrums aus (42.4) als Lösung der TSCHEBYSCHEFF-Aufgabe 
mit der Referenzabweichung (42.2) als gesuchtem Minimum. 


3. Schritt. Ermittlung einer aus der alten Referenz ausscheidenden Zeile Ay: 
 Aufsuchen einer speziellen Lösung des Gleichungssystems 


a, = m: U* (m; = (km Un. -Mrta)) > =. - . . (42,5), 
‚ Bestimmung eines Index /, für den entweder gleichzeitig 
| Ah=0 und — sn mh DDr une . a 
' gilt oder — falls ein solcher Index nicht existiert — für den das 
: A ki,k 
max (—T- sgn nk a In ht a Te a 
(1:2, #0) Ai 


angenommen wird. 


4.6. Die an Formel (42.6) anknüpfende Ausscheidungsvorschrift bedarf freilich einer 
ı Ergänzung, um „Referenzenzyklen‘“ zu vermeiden. Wir wollen von einer „ernsthaften 
ı Entartung‘“ sprechen, wenn es mehrere Indizes i mit der Eigenschaft 

| kV und — sgn nr Pie > 0 

gibt. Falls die unmittelbar vorangehende Referenz mit keiner ernsthaften Entartung verknüpft 
ist, falls es bei dieser höchstens einen Index oder gar keine Indizes mit jener Eigenschaft gibt 
wollen wir von einer „lokal erstmaligen ernsthaften Entartung‘ reden. Eine lokal erstmalige 


ernsthafte Entartung könnte der Anfang eines Referenzenzyklus sein, wenn man unter den 
Indizes gleicher Eigenschaft immer je einen willkürlich auswählt, 


L. BITTNeR, Das Austauschverfahren der linearen TSCHEBYSCHEFF-Approximation ... 249 


Die Matrix mit einer solchen Entartung möge 


Aw "e-- 0 
Fi = 
* An, er Te . . . . . 2 
W=|1: und die Vorzeichen ihrer Zeilen mögen 
- : 
un 0 
Ao,ıı Besi 
heißen. Wir ergänzen W,, a... , @r+ı In irgendeiner Weise durch @,+2,...., @, zu einer 
Permutation von 1,2,...,n und ordnen den Zeilen Co u a A, irgendwelche Vor- . 
zeichen ß}+3,...fn willkürlich zu. e sei eine hinreichend kleine positive Zahl. Wir defi- 
nieren damit ein zu Problem 4.1 gehöriges e-Problem: 


3(e)-q = Maxl 
für 7-30 + LH 
j-1 


I 


3 n+1)to; 


a + e2r+1.(0,,1), ME benrı. (43). 
ac I (n+1)+0; 


2 


Man findet durch die gleichen Umformungen, wie wir sie in 4.2 und 4.3 durchgeführt 
haben, daß eine zulässige Lösung 3(e) des e-Problems eine Lösung des Systems 


. I E n i I 5 antl 
3(e)- m\=LEß Ar 2 et) u; > een =&eE . . (44) 
—g j-1 j=1 j=0 
ist?); insbesondere ist eine zulässige Basislösung 3*(e) die Lösung eines Systems der Gestalt 
hr i nn; % - £ 2n+l 
3*(e)- m) = Ra, + Lerti— Pi) a, (ehe (A); 
er j-1 ; i= 
A\* 
0, | bezeichnet hierbei also irgendeine beliebige auf W5 folgende oder mit W$ übereinstim- 
—Wy 


mende Referenz des &e-Problems und ist deshalb auch eine Referenz des Problems 4.1. 
Es sei %* die dazugehörige Basis und es werde 


EN bes ae IR a Eee El: 5), 
gesetzt. Man findet aus (43), daß 


5 n 
+Letih, n + eintlzrr . .. (46) 
n+1)+0; 


n 
*(e)— 3’ NV Ep 0 
P 1—P; : 
’i—ı ar D+e, j=1 


2 


ist. Nach Substitution dieses Ausdrucks in Formel (44°) und darauffolgendem Koeffizienten- 
vergleich in der Unbestimmten & erhalten wir endlich füri=1,2,...,n 


A\* A\* 
d=b; De vs al und — = %ııti Un.|: > Dnzıdi® er =1 (47). 
—HY —U 


Wir nehmen nun mit den Komponenten der Vektoren d;, Yu+ı+i, 5*(e) und 3* die nämliche 
Umbenennung vor wie am Anfang des Abschnitts 4.5 und können dann beispielsweise schreiben 


Pı 
EEE VE a: vier =1. 


Pr+1 


3) Das System von Nebenbedingungen (43) läßt sich nämlich mit Hilfe einer gewissen Permutations- 
matrix ® in der Form 3’()-3 = (0,1) +E ®% darstellen. Daraus ergibt sich: Ge) -ERS = (rl), 
d.h., der Vektor 3’(e) — € ® muß den im 1. Ergebnis formulierten Bedingungen genügen. 
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| Mit Hilfe einer speziellen Lösung m; = (fi +++» r+1,,)‘ des Gleichungssystems 
| a) = mi-W* EL 


| kann man also den Vektor dv, und Dn+i+i in der Form 


ßı Pı £ 
1 mtl,  ne| .-m+sl) .. (4) 


Pr+ı Pr+ı 


v1, = 


. 


darstellen; die Skalare s# und s; lassen sich aus der Normierungsvorschrift b; - &+1=1 bzw. 
Dazızs" &r+ı = 1 errechnen. 
Statt (46) schreiben wir mit Benutzung der Formel (49) auch 


Pı ri] hr. g° N 4 —B} 2n+ı r-| 
o=| aa tele ++ Kerl mu +5 Er = 
r-+1/ (5 )- 


Die Serie aufeinanderfolgender Referenzen mit Entartungen, die wir mit Hilfe des e-Pro- 
blems auseinander herleiten, beginnt mit der Referenzenmatrix U5 und endet mit der Referenz, 
für welche die Indexmenge {i :A; = 0, — sgn nx Pi wi,e > 0} (lokal erstmalig) leer wird. Bei 
einer jeden Referenz dieser Serie — mit Ausnahme der letzten — haben wir daher zur Ermittlung 


Fe SER ZEN EN . 222Me) 
der auszuscheidenden Zeile a,, das min bzw. min Fan 
(1:9,5>0 Dir Pin+irk 
festzustellen, für den das min angenommen wird. Da es Indizes mit: z> = IA:| —0( und 
— sgn x Pit,» > 0 gibt, ist von vornherein klar, daß das gesuchte min nur für einen Index 


iel={i:4=0, — sgn ng + Bi Wi, > 0} angenommen wird. 


zu bilden und den Index I 


Wir haben uns also zu fragen, für welchen Index i =! das 


2 ar ” 
1 0 +5 A8 ;- 
min —— p.|& EP} ui,a; — DE" By hie, 
viel D;,k . j=1 ji=1l 
bzw. — in) = ee 
: zE iel — sgn nr Pi Mi,x 
min ———— (51) 
ieI Di,n+1+k 
'n n 
5 Ki,o; SR Ko; 
— — max (sen ne |, EPß}- I — Di) entigg. 3 
ieI j-l ik j=1 Kir 


angenommen wird. Dieser Index ist eindeutig bestimmt und kann offenbar in der folgenden 


. r Li 
Weise ermittelt werden. Man ermittelt die Indizesie I, für welche das max (son n«BR- Er > ) 
| ie i,k 
angenommen wird; diese Indexmenge heiße I,. Besteht /, aus mehr als einem Index, 


ik 


‚ so stellt man die Indizes fest, für welche das max (son n« P2- a) angenommen wird; 
ieI, 
| diese Indexmenge heiße I,. So fortfahrend, gelangt man schließlich zu genau einem Index I. 


5. Tabellenmäßige Fassung des Austauschverfahrens 


Bevor wir an einem Beispiel die Handhabung dieser Ausscheidungsregeln demonstrieren, 
wollen wir noch auf die Verwendung der von STIEFEL eingeführten reduzierten Tableaus hin- 
weisen. Es seien unter den Zeilen RR ET nr einer Referenz z. B. gerade die r ersten linear 
unabhängig. Jede weitere Zeile a; kann dann eindeutig als Linearkombination 


= 50, 7 dar Ant hr ti ee KR ae ee 
insbesondere 
RR nn: A, a, + Aa 0, >» /r An, 
dargestellt werden. Die Größen Kij ergänzt um 41,5 = 0, und A, ergänzt um ,ıı = —1, 


können als spezielle Lösungen von (48) und (36) angesehen werden. Aus (52) folgt: 


N (& 2 Mij cu) + z MiiNv; > Re — (or: HR e) +24 N;: 
se den = i-1 


i=1 


H 


a 


me A x u EEE FT W Pe #7 
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\ 


Diese linearen Abhängigkeiten notiert man zweckmäßiger Weise in Form einer Tabelle 


i . | nv; Bi ... N9,+1 > ale 9 n= ER - NZ 
nn nn nn nn 
1 Mm Pı A Mi,j Hı,k 
I Nu; Bı hi u; 1,k 
% N Pr Ar Ur,s Ur,k 
or j ae == 
ar Nor 41 Pr+1 En —1 an: 0 Je 0 | 
r+1 14 ‚ Lu \ 
— 24 C; G—% M,joo, Ck — 2 Mi,k %; 
i=1 si i=1 
— Ir +1 => = 9% 


In Gedanken möge man sich diese Tabelle um r weitere Spalten mit den Überschriften 
N No +++, ergänzt. denken, in denen die Koeffizienten 11, = 1, u2,,=0,.:..:,Wı4, = 0; 
kun, =, un, = 1 us, =0,...,Hr+ı,., = 0; der linearen Abhängigkeit der Zeilen a,,0,,.-- 


von diesen Zeilen a,,4,,»..-; a, a4 selbst notiert stehen. 
Aus einer solchen Tabelle liest man z. B. sofort = sgng,,„h= mal ab. Die 
1+2]4| 
i=1 
Funktionswerte der. n; für die Referenzzentren lassen sich wegen n,, = — ß;- h ohne viel Mühe 


errechnen. Von großem Vorteil ist der Umstand, daß man zur Herstellung der darauffolgenden 
Tabelle nicht Gleichungssysteme der Gestalt (52) auflösen muß, sondern lediglich eine JORDAN- 
Elimination auszuführen braucht; wenn etwa gar 6, — 0 vg die auszutauschende Zeile ist, 
erübrigt sich auch diese Elimination. Sonst muß man eben n; gegen n,, (1 S!= r) austauschen 
— vorausgesetzt, das „pivot“ u, ist #0. u. ist gewiß + 0, wenn die Indexmenge I nicht 
leer ist; 1 „ könnte also nur dann verschwinden, wenn die Ermittlung des Index ! nach Formel 
(42.7) erfolgen muß. Das max der Formel (42.7) ist infolgedessen selbst = 0 und wird wegen 
A = —17#0 und 441,2 = 0 offensichtlich auch für den Index i=r-+ 1 angenommen. 
Um Rechenarbeit zu vermeiden, wird man daher ER “ als auszuscheidende Zeile 0, betrachten. 


Über Jorpan-Elimination informiert ausführlich [2]?). 


6. Ein Beispiel 


Gegeben ist das System von linearen Funktionen 


Muh 3. +4 +10, 9=%at+% +2, —7, 
13, = +12, —5/2%, — 2, +% +32, + 37/2, Mt —%+6, 
Not —u+3; 1 = — 2. —5/2%+ 13/22, + — 2 —920— 15/2, 


17=2% +%9 7%, +2, —22, +28, NY =. +52, — 1222, —2, +52, +32, + 11/2, 
mr Ihr th: 


Der Rang der homogenen Teile dieser linearen Funktionen ist r—=4. Die Zeilen a1, a3, 
a4, a, sind linear unabhängig und bilden zusammen mit ag eine erste Referenz. Die folgenden 
Tableaus enthalten die Ergebnisse der schrittweisen Rechnung. Die mit „f. RZ. beschriftete 
Zeile und Spalte dieser Tableaus gibt die Werte der Funktionen n,; in einem Referenzzentrum an. 
Die Referenzabweichung der ersten Referenz ist 0; wir setzen hier willkürlich 7 = + 1 fest 


4) Diese Eliminationen sind nicht auf die Zeilen 7,, ,, anzuwenden. 
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Bw 


TI »+, ml Peer 
Il er HER Nee u His 22 na 
Mi,s\ en max (resume 9) =0 für i=1=5, u= 
max (a sgnna 2) =1 ri (4;#0) 5 
a i 
| aa | M | | min m mERZ. 
mn || mm mim) m Im m RZ mm 
Rh | | he 
pr 1 —ı 0 21 —1 i 
LE LE a 5 BE BE 1 alt ae o|-ıl ı1-ı2) ol —ı 
ml len da ılooa em lee Iron 
3 | ns ww u N HE ee Fe ee 
\m |#+| 12) 2) 940-1 |olo!-ı 4| » |-]| 4 -—ı| 0] 0] oje, ı, 
= Ar 4 Te a es Er EN ae 3 Eee TE a 
‚5 1 az Wa 0/10 0 0|0,0 L 5 ns = 0| ol—1] 170, - 
l - 
Be neh Sa 
! äyt I a ! 
ERZ.B-.) 31 0-39", me £.RZz. | ıla=ı| 0] o| 2] m 
T=4+1, |n|l>h, P=—sgenn =—l1, z=+1, ml>h, B=—-sgn=—1 
1=41,2, 38; I={1,33, L={l=2, n=2 
Anfang des e-Problems mit &, =1,...,@%=9, 
Pr. cut = 1, a lie Ba— eh — LT: 
De A 
| Be 
im; Alm mim|ı m.| m nn fRZUOGE m || Mm nm | m| m |ERZ 
EEE. 91-1 ı|m +) 4/5|—6/5| 12/5| 1] 14/6] —ı 
3 16 IP EEE ee Be et Meng Bi = 
& N ea 0 1 2. un 2 = | 0 | °| 0 —ı| 0 Ba, 
EL ENT 3 | —I-35| 25l—a5| 0|—8%| 1 
| m |-—|—| ıl 0o| 2|-—3| 2 4m |-j—15| Mala) oje ı 
5| m» |-|—2| 2| 0| 2|—| 7 5|m |-]-265|—275| 45] 02%) ı 
| Bin nA 9 | | 951—6/5|+7/5| 2| 296 
£.RZ 4 ld E Inn | Rz. |-ı5| 461-875 k=1l—ıjs| vi. 
r=+l, |n|>h, B=—sgsann = —|1 Tesrl, In >h, =—gnn =+1 
rel) = I=%4,. Lech. n=8 
— — - 
im Bl Mm ne | m \mi 7 |ERZ. 
R ’ ER Ende der Rechnung, da alle |n5| < Ah sind. 
098) re] BT 2—5/3| —1  Diezugehörigen Referenzzentren, Lösungen der TscHE- 
a I a nn 2 1 2 # | 1 »ysomker-Aufgabe, haben die Koordinaten 
Ba 22 ze Ba VER BE 
2» || 0 | 0 |-ı[ Jo IT 2 = — 134/38 — 7/3 25 + 8/8 20, 
4 N6 nn > —1/3 —2/3 0 0 —4/3 +1 a = 14 — % — Amar, 
el UN 6 ni) « 
na BRLOF 212 1 WB] ve u Se 
> 4/3 2/3| 2 8 7/3 | u=— WW +13: 0 +4/3 2, 
{.R} ) 0 |A=ıl U 0 |vı. %, 2; = willkürlich 


—_ ww ng - 
g % IR 5 2 F - 
4 ER ; j r & " J 
v . 
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7. Die Gesamtheit aller Lösungen des Tschebyscheff-Problems 


7. Bei nicht-erfüllter HAArscher Bedingung ist die Lösung der gestellten Aufgabe bekannt- 
lich nicht mehr eindeutig. Um unter allen möglichen Lösungen eine solche, die noch gewisse 
andere Vorzugseigenschaften besitzt, auswählen zu können, braucht man eine Darstellung der 
Gesamtheit aller Lösungen. Diese Gesamtheit ist identisch mit der Gesamtheit der zulässigen 
Lösungen x des Systems linearer Ungleichungen: 


det, — n.() sh a ER (53); 
h bedeutet hierbei das (ermittelte) TSCHEBYSCHEFF-Minimum. Durch die Ungleichungen (53) 
wird bekanntlich ein konvexes Teilgebiet des n-dimensionalen Raumes, ein Polyeder, definiert. 
Im Falle r = n ist dieses Polyeder beschränkt, im Falle r <n dagegen nicht. Um die Unter- 
suchung auf den ersten Fall zurückzuführen, führen wir mit Hilfe von 


De en Ze Mi arten Dar EHER (54) 


neue Koordinaten ein — ähnlich wie in (6). v, v3... ,», seien gewisse aus der Reihe 1,...,n 
stammende Indizes mit der Eigenschaft, daß R(a,, ,...,q,)=r ist. An die Stelle der Un- 
gleichungen (53) treten nun Ungleichungen der Form 


HER, —L;4) sh a ee De ee! (53%): 


Im Raum der Variablen (y,,y,...,y,)’ =) ist das Polyeder der zulässigen Lösungen 
von (53’) beschränkt, denn das zu (53°) gehörige System homogener Ungleichungen 


:Yy)=0, 6) = 0 (=12,...,n) 

enthält u.a. die Ungleichungen y, <0,...,.y, <0,—yı S0,...,—y,<0 und kann somit 
nur die Lösung yy=y=---=y,—=0 haben. Nach einem Satz über lineare Ungleichungen 
(siehe [5], Theorem 1 u. Corollar 1 B) ist die Lösungsmannigfaltigkeit von (53°) beschränkt u. 
kann dargestellt werden als konvexe Hülle ihrer „Ecken“. Eine „Ecke“ ist eine spezielle zu- 
lässige Lösung } von (53°), für die bei wenigstens r Funktionen mit homogenen Teilen vom 
Range r die strenge Gleichheit (= statt <) gilt. Die letzte Referenz des Austauschverfahrens 
liefert solch eine Ecke. Kennt man also sämtliche Ecken, so kennt man damit auch alle mögli- 
chen Lösungen von (53°). Eine naheliegende, aber sehr umständliche Methode, sämtliche Ecken 
zu ermitteln, besteht darin, daß man alle möglichen, aus r Gleichungen vom Range r bestehende 
Teilgleichungssysteme von (53) bildet und nachprüft, ob die zugehörigen Lösungen den übrigen 
Ungleichungen genügen. Man erspart sich aber gewiß viel nutzlose Arbeit, wenn man zu jeder 
Ecke zunächst die „benachbarten“ Ecken aufsucht, die mit jener unmittelbar durch „Kanten“ 
verbunden sind. Zur Konkretisierung dieses Gedankens nehmen wir an, daß 9® eine gewisse 
Ecke und : 


Kar &,®) =h, Ye, ... Yp' Sy) = h pzn, (55) 
y;-5®) <h (für alle übrigen Funktionen) ’ ’ 
das zugehörige Ungleichungssystem (53°) sei; y1, - . - , Y), y; bedeuten Vorzeichenfaktoren. Wegen 
y;:&(y®) < h gilt auch noch y; -£;(y) < h für y — y| < 1. Daraus folgt, daß die hinreichend 
nahe bei 9) gelegenen zulässigen Lösungen )) von (53°) allein durch die p Ungleichungen 
yı'&W) sh, 2 erh EEE . (56) 
bestimmt werden. Aus (56) und (55) folgt überdies 


y 090) = 0, HMI HN EO .. 6. 
Offenbar kann jede Lösung von (53°) mit Hilfe einer Lösung y* von (56) folgendermaßen dar- 
gestellt werden: y = y® + 1(y* — 99), O <t< t(y*). 

Hat also (56) keine weiteren zulässigen Lösungen mehr, so gibt es auch keine weiteren 
für (53°). 

Durch die Ungleichungen (57) wird ein „Kegel‘‘ mit der Spitze in 9!) festgelegt. Es ist 
nicht ausgeschlossen, daß sich unter diesen Ungleichungen einige überflüssige Bedingungen 
befinden, Bedingungen also, die immer erfüllt sind, wenn die übrigen Bedingungen gelten. Wir 
wollen annehmen, daß die Beziehungen 


Ya+tı' Cu le) = Nm. Yp’ a 6) — y®) <(0 (q =, p) 
immer erfüllt sind, sobald nur 
rend) EFERAGH 


gilt, daß diese q Ungleichungen aber nötig sind, um denselben Kegel wie in (57) zu definieren. 
Mit 5b, werde der Normalenvektor der Gleichung £;(y) = 0 bezeichnet. Nach einem Satz von 


N 


._ 
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Farkas (siehe [6], Theorem 3) ist jeder der Normalenvektoren Yy+ı Bg+1 «+» » %» b, positiv 
linear mit ‘Hilfe der Normalenvektoren y} bj, - - +» Yg a ausdrückbar, liegt also in dem von 
y161-- 946, aufgespannten Kegel. Daß auch umgekehrt eine lineare Ungleichung aus (57) 
immer mit den übrigen Ungleichungen erfüllt ist, wenn sie positiv linear von diesen abhängt, 
ist unmittelbar zu sehen. Zur Aussonderung überflüssiger Ungleichungen hat man also die 
positiv linear abhängigen Normalenvektoren ausfindig zu machen. 


Weil die Vektoren Yy4ı bat ++. ,7%pdp von Yıb1...,Ygbg abhängig sind, folgt, daß 


schon 

RyıbYabo.:-,Yb)=T also fe En En 1 > (59) 
ist. Bekanntlich läßt sich mit Rücksicht auf (59) aus der Existenz einer linearen Beziehung 

ybehydbtie pt + la Ya da 

mit nicht-negativen Koeffizienten f,, 1, ... f, immer eine lineare Beziehung y b = U Yo, da, + 
BY tr EHNEH..- SG < q mit ebenfalls nicht-negativen Koeffizienten 
t,Bs...,t, und linear unabhängigen Vektoren y,, ba, - - - »Yo, Do, herleiten. Daraus darf 
man schließen, daß jeder der Normalenvektoren überflüssiger Ungleichungen positiv linear 
mit Hilfe von je rlinear unabhängigen Normalenvektoren der übrigen Ungleichungen aus- 
gedrückt werden kann. | 


Der Einfachheit halber sei b,, b,, ... . , b, ein solches System linear unabhängiger Normalen- 
vektoren. Mit Benutzung der Tabellen der Nr. 6 erhält man die linearen Abhängigkeiten 


yzkıM FkiaNnt + er th: r<i5Dp ec te ee (60) 
und damit 
ybı = Yyıyılliı Yıdı + YıyalliaYabat "+ Yiyrlür’ rd - - » - (60). 


Man darf mithin alle diejenigen Funktionen außer Betracht lassen, für die 


Yyıliız 0; Yyıkiazls 1.» Yyrbunl nn kb 


gilt5). Danach werde eine andere Größe n;(r <i <p) mit Hilfe einer JorpAan-Elimination zu 
einer darstellenden Größe gemacht und wiederum die Aussonderung der überflüssigen Funk- 
tionen analog zu (61) vorgenommen; usw. Auf diese Weise können alle überflüssigen Unglei- 
chungen des Systems (57) erkannt und entfernt werden. 


Sämtliche „benachbarten“ Ecken des Polyeders (53’) liegen auf den „Kanten“ des Kegels 
(58). Eine Kante von (58) ist eine eindimensionale Lösungsmannigfaltigkeit, für welche strenge 
Gleichheit genau bei einem Teilsystem vom Range r—1 gilt. Daß es solche Kanten gibt, 
wofern (58) außer für 9 auch noch für andere ) Bestand haben soll, ist klar, denn die Gesamt- 
heit aller zulässigen Lösungen 9 —9y® von (58) kann als konvexe Hülle eben dieser Kanten 
dargestellt werden (siehe [6], Corollar 1 B). Zur Ermittlung sämtlicher Kanten von (58) können 
auch wieder die Tabellen der Nr.6 verwendet werden. Es sei nämlich z. B. der Rang der ersten 
r— 1 Funktionen des Systems (58) gerade r— 1. Die allgemeine Lösung des Gleichungssystems 


na) =0, EM), Yard) = 0 
läßt sich mit Hilfe einer speziellen Lösung y* — )® in der Form y — y® = 1(y* — 9) angeben. 
Unter den q Funktionen von (58) gibt es wegen (59) sicher eine Funktion, welche zusammen 


mit den r— 1 voranstehenden ein System vom Range r bildet; es sei y,£,(d — 9) eine solche 
Funktion. 


_ Offenbar ist es durch geeignete Vorzeichenwahl stets möglich, zu erreichen, daß 
Yrörky* —Y) <0 gilt. Die übrigen Ungleichungen von (58) verlangen 


vs) <S0 w<isg, 


wenn )* — 9) zulässige Lösung sein soll. Aus (60) ergibt sich 


ar) = 2 yıyania ya d* YO) = Ya Yrhinr yet — HR). 
Soll dieses < 0 sein, so muß y;y, u, Z0 i=r 

< s IYyyrkir zz (i=r+l1,...,g) gelten; trifft das zu, so stellt der 
Halbstrahl 9 = y® + £(y* — 9®) (t > 0) eine Kante von (58) dar. 


5) Aus der Beziel 30’ } ; Su 
) ziehung (60') kann man auch gleich noch erkennen, ob etwa ein Vektor Ybbs lSHSr) 


positiv linear von y;b; und allen anderen Ybs 1SsSı ä i 
sb 1&Esszr,s=+ s,) abhängt. Das ist nu 
außer »s; #i,s, alle anderen 5: 1, < 0 sind. » 0) g £ nur dann der Fall, wenn 


ei u BE u» ei Ye Fe NE 0 En Ze 
r - 


Be 
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— 


Offenbar kann man also alle möglichen von 9 ausgehenden Kanten dadurch erkennen, 
daß man von allen möglichen Funktionssystemen Na ++» No, (0, = gq) des Ranges r ausgeht, die 


zugehörigen linearen Abhängigkeiten der Form 


m = 2 Mies 1% RE BE (62) 
betrachtet und feststellt, ob bzw. für welchen Index s = ! die Zahlengrößen 


Yıyaz Wi 0 VE ee ea Ark (63) 


sind; falls (63) zutrifft, geht von 9® eine in den Ebenen y,, - er Y—9®9)=0 (s =# I) gelegene 
Kante weg?). | 


Nachdem man r — 1 Ebenen gefunden hat, in denen eine Kante verläuft, ist es nicht mehr 
schwierig, die Ebene y; - 3,(y) = h zu entdecken, welche die Kante in einer benachbarten Ecke 
des Polyeders (53°) schneidet. Der einfacheren Schreibweise wegen nehmen wir wieder an, daß 
die Kante in den Ebenen der (r— 1) ersten Funktionen von (58) liege. Von dem Halbstrahl 
y=YyV +1(y* —9®) (> 0) gehört nur dasjenige Stück zum Polyeder (53), für das auch 
die Ungleichungen 


3-0) = yo) Hy ON) Eh,  jeti:y.G0)<h 
erfüllt sind. Da y;-£;9* —y®) = y5 14,,&.(y* —Y®) und nach voraufgehender Vereinbarung 
6. — m) <O ist, erkennt man, daß eine solche Ungleichung im Falle y; y, u,,, Z 0 keine 
weiteren Beschränkungen mit sich bringt; aber im Falle y; y, t;,, < 0 muß 
 hon:509) 
Yin &r* Hy) 
gewählt werden. Es tritt in der Reihe der betrachteten Ungleichungen y; - £;(y) < h wenigstens 
ein Gleichheitszeichen ein, wenn £ gleich dem min dieser Schranken genommen wird. Man 


findet daher eine die Kante schneidende, mit dem benachbarten Eckpunkt inzidente Ebene 
durch Ermittlung eines Index / = k, für den das 


h— y3: 09) 1 un EL E 


min : Z en le max - 
g: v; Ihr Yilir 6 (y Ir y®) Yr & (y* se N) dr Yyrli,r = YiYrkir 


also das max der rechten Seite angenommen wird. Die Ebenen y,&.() =" Y—-ıd&r-ı()) = 

y1&4(9) = h bestimmen den gesuchten Eckpunkt. Seine Koordinater y, = 1, Yg = Ma ++ Hr = N, 
kann man einfach dadurch erhalten, daß man mittels einer JorDAn-Elimination n, gegen nı 

austauscht und die zum =Yyıh,...,9m-ı = yı—ıh, nr: = yrh gehörigen Funktionswerte von 

ER errechnet. Auf ähnliche Weise hat man vorzugehen, um alle zu 9) benachbarten 

Ecken zu bekommen; von jeder dieser erhaltenen Ecken ausgehend, kann man ebenso die dazu- 

gehörigen benachbarten Ecken aufsuchen; usw. Wichtig ist, daß nur die Tabellenwerte «; , und 

zu deren Herstellung JOoRDAN-Eliminationen vonnöten sind. 


V,9y2,...y® mögen das System sämtlicher Ecken bilden. Jede zulässige Lösung ) 
von (53) ist darstellbar als 


y=arıyP Hyd + ...299, Zw;=1, u EN ee lo 
=1 


Eine partikuläre Lösung r des Gleichungssystems 


(7,8); N,(2); Hide; 7,8) =E yo 


heiße r®, und die allgemeine Lösung x des homogenen Systems 7,,(£) = + ,(0) = 0 heiße tr. 
Die Gesamtheit aller zulässigen Lösungen des Systems (53), d. h. die Gesamtheit aller Lösungen 
des TSCHEBYSCHEFF-Problems, läßt sich nun in der Form 


gend +... 0 + LE, Del, Sl, une» (69) 
i=1 


darstellen. 


6) Die Koeffizienten w;,s in (62) sind selbstverständlich von der Kombination o,, 03, ..., or abhängig. 


Beispiel: Aus der letzten Tabelle der Nr. 6 entnimmt man, daß ww 3 

W=W=-P-—1, P-P--T Won 
eine Ecke und das System . N 
Sl, -ysl, „Sl, »Sı, Ban+2sı, m-2m 22m +31 
den Kegel (56) bzw. (57) bildet. Die zweite Ungleichung ist entbehrlich). Ein Austausch yon 3 
4 gegen n, führt auf die Tabelle: «4 


al na] ze) Se] 8} 
u 


Unter den übrigbleibenden 5 Ungleichungen gibt es also — wie die Spalte für n, (siehe®)) lehrt — 
keine überflüssigen mehr. Mit Hilfe der Formel (63) (hier i=5, ,=3,,=-7,,=a=6%, 
04 = 9) erkennt man, daß in den Ebenen —y =1,1,=1, 7% = 1 eine Kante verläuft. Im 
Anschluß an Formel (64) hat man zur Ermittlung der 4. Ebene unter allen, nicht zu den obigen 
6 Ungleichungen gehörigen Funktionen: +7, EN» +2» + Ne — Ns, — Ns — N» E Ns — Ne 
zunächst diejenigen aufzusuchen, für die eine der Relation y; y, #,.< 0 entsprechende Relation 
gilt; (die «;, sind der 3. Zeile der vorigen Tabelle zu entnehmen.) Das sind die Funktionen 
+» +» —Ne +Ns Das max (— 3; — 3/2; — 2; —3/4;) = — 3/4 und wird für +7 
angenommen. Der Austausch von n, gegen nz liefert die Tabelle: 


. 


BIENIEITSKERT 


m|-al-ıel ı | ı auLn.D 
m | 12 0 | —ı2|! -ı2 —32° + 
I 1a} al. 0 7 ee 
ns 0 0 I —ı I—ı I—2 | + 
Isa! -ı2]) zel 3e| as | 
Die 2. Ecke lautet daher y9— (u =—1l, ®@=—1, A= a Sn un 1) 2 


Jede andere der 3 weiteren Möglichkeiten führt auf dieselbe Ecke. 
Die Gesamtheit aller Lösungen des Beispiels läßt sich infolgedessen durch 


— 134/3 — SPN — BB, +8Bx 
14 29/2 —. —n 
a u = ng 0<x<siı, 
Pe g 1 Br n ES 5 6 b 
! — 1/3 +( a, — 1/31] 13%, +43 x; x, willkürlich , 
0 Is 
\ 0 0 Te 
darstellen. 
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ZAMM 41 (1961) Heft 6, Seite 257—261 
Der profil-ebene Spannungszustand *) 
Von WOLFGANG VOckE | 


An dünnen Scheiben variabler Dicke wird ein profil-ebener Spannungszustand definiert. Dieser 
hängt im allgemeinen von der Poisson’schen Konstanten ab. Aus den Sonderfällen, bei denen Unabhängig- 
keit von 1/m vorhanden ist, resultiert eine ‚‚verkürzte‘‘ Differentialgleichung von einfacher Form. Bei ent- 
sprechend begrenzter Profilbreite stellt diese eine gute Näherungsgleichung für beliebige Profile dar. 


For thin disks of variable thickness a profile plane distribution of strain is defined which in general 
depends on the Poisson ratio. In the special cases (independence from 1/m) a ““shortened’ differential 
equation of simple form is obtained which is a good approximation for arbitrary profile provided the breath 
of the profile is suitably limited. 


Ha TOHKUX IIIaCTUHKax IIepeMeHHOÄ TOJINHHBI OUPeNeNAeTcH „„HPOBUAIBHO-INOCKoe“ 
HAIPASKEHHOE COCTOAHME. B oÖMeM cıLyyae OHO 3ABUCHT OT IWOCTOAHHOU IIyaccona. Jun 
CHEeIMAJIBHEIX CIIyYaeB, Te HET 3ABUCHMOCTH OT 1/m, NudhhepenmmanbHoe ypaBHeHnue Öe3 
ıpa6oüi yacTu UMeer IpocToü Bun. IIpu COOTBCTeTBeHHO OTPAHHYeHHOH IUUPUHE ITPohuma OHO 
ABJIHETCH XOPOIHUM IIPHÖJHFKEHNEM YPaBHeHUA MIA JIOOBIX IIpodhnseii. 


Definition und Differentialgleichung 


Eine Scheibe mit variabler Dicke h werde symmetrisch zu ihrer Mittelfläche z = 0 bean- 
sprucht (Bild 1). Die Dicke h wird gleich a - f(x, y) gesetzt, wobei a einen konstanten Faktor 
und /(z, y) die Profilfunktion darstellen. 5 
Da die Scheibe beliebig dünn sein soll, ver- 
schwindet bei der oben vorausgesetzten sym- 
metrischen Beanspruchung die Schubkompo- 
nente in z-Richtung. Weiterhin soll an allen 
Stellen der Winkel zwischen der z-Achse und 
der Normalen auf der Profiloberfläche klein 
sein. Den so definierten Spannungszustand 
wollen wir als profil-eben bezeichnen. 


Die Gleichgewichtsbedingungen 


(% 


N 


| 9% 2 
RER a Bild 1 
(ey 3 Ds y sr (Try z Dsz —=0 
lassen sich durch eine Spannungsfunktion erfüllen. Setzt man nämlich 
F, yn 
= 
Ser 
ö F, 
en VE a 06 1 N En ee ES) 
und Ba 
Try = tr 


so sind die Gleichungen (1) erfüllt. Die Kompatibilitätsbedingung liefert unter Beachtung des 
Hooxsschen Gesetzes die Differentialgleichung für F: 


ee 


bzw. mit Einführung des A-Operators in der Ebene: 


AAF— KANNE E fr - (AF), 442 
. 1\ N 
a [E22 Rn (4). 


they 2 far BF 2tf- hey 2 hart Fızu) =0 


*) Erster Teil der von der Fak. f. Maschinenwesen der T. H. Dresden angenommen Habilitationsschrift. 
(Referenten: Prof. Dr. A. WEIGAND und. Prof. Dr. H. NEUBER) 
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Im Grenzfall f— konst. bekommen wir den ebenen Spannungszustand mit AAF=O. 
Die Arkysche Spannungsfunktion ist also auch ein Sonderfall der profil-ebenen Spannungs- 
funktion. Weiterhin sieht man, daß die profil-ebene Spannungsfunktion von der Poıssoxschen 
Konstanten im allgemeinen nicht unabhängig ist. 


Die Randbedingungen 


Wenn wir mit X und Y die Komponenten der resultierenden Randspannung bezeichnen, 
so liest man aus Bild 2 die Randbedingungen ab 


Y-F.ds anne: 6) 


Iy 12 of de—ny-f-dy—Y-f-ds=0 


N bzw. nach Integration über einen Randstreifen 
X:f.d5 me di er Ey und unter Benutzung der Gleichungen (2) 
£ y es \y 
| Zu Fa. PR =—(R ht ER (6) 
1 P,=(F, a — (F,2}ı 


| j 0,-F.dx 
= Die Formulierung der Randbedingungen 


Bild 2 bleibt dieselbe wie beim ebenen Spannungszustand. 


Sonderfälle der Profilfunktion 
Das Keil-Profil 


Falls f nur von y abhängt, also ein beliebiges Keilprofil (Bild 3) beschreibt, vereinfacht 
sich Gleichung (4) zu 


8 2» Pu —hewf) np 
MIR a (Pr n )=0 zur RR 
h=0-fir) 5 
GEEEEEEGE EGE EEEEEEEE 5 
SRRIIRIIIÜN . 
N 
N N 
Ba Y 
Bild 3 .- 
Y 


Das rotationssymmetrische Profil 
Wenn wir Polarkoordinaten 
E=T.Cosp, x 


y=r-sinp 


einführen und / nur von r abhäne \ i RK 
abhängen soll, also ein beliebioes ne 2 . 
i R £ als biges Rotationsprofil z-A 
s Br Pe en ZN f g ße um die z-/ a 
schreibt (Bild 4), erscheint Gleichung (4) in der Form: p 15 chse be 


: 1 
MAP 2 BE ae 
(AF) j (AF),, 7 

ef 


+ Uhn—2. mo) Seel] =0 


lea hl Enten 


m“rT 


8) 
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Profil-ebene Spannungszustände, die von 1/m unabhängig sind 


Will man diejenigen Profile finden, für welche der profil-ebene Spannungszustand bei 
jeder Belastung unabhängig von der Poıssonschen Konstanten ist, so müssen in Gleichung (4) 
die Bedingungen 

Elke =0, 


ER a EN NE In (9) 
le Zei =0 
erfüllt sein. Das Gleichungssystem ist in x und y symmetrisch, und die Lösung lautet 


ee Ar en) 


Neben dem ebenen Spannungszustand (f = 1) ist also auch der Spannungszustand am 
hyperbolischen Keilprofil (Gl. (10)) von 1/m unabhängig. 

Umgekehrt bekommen wir bei beliebigem Keilprofil einzelne Spannungszustände, die 
ebenfalls von der Poıssoxschen Konstanten unabhängig sind. Setzen wir in Gl. (7) 


1 0x Te 0 D 
dann können die Lösungen nur folgende Form haben: 
4 8 
zu) er 2 Nu)e 2 REF EEE (11): 
E = u=5 
Sie lauten im einzelnen (nach Belastungsfällen geordnet) 
Fa =K, +R.y (a) 
en R,:.It Kr: us] (DIE Te A ee (12). 
F, = Kr + Kr r-y (e) 
Fa, Km IH Remy l (d) 


Für Zug, Druck, reine Biegung (12b), Schub (12c) und Biegung mit Querkraft (12d) ist 
die exakte Lösung identisch mit der elementaren Formulierung, wenn / = /(y) und o, = 0 gesetzt 
werden (Bild 5). 


Bild 5 


Schließlich lassen sich bei zwei bestimmten Belastungen die Bestimmungsgleichungen für 
f angeben, so daß F unabhängig von 1/m ist. Die Profilfunktion / hängt dabei von x und y ab. 
Im ersten Fall setzen wir 


F=2(@+9) BEE EA ER E NN 


Das entspricht einer allseitig gleich großen Normalkraft/Längeneinheit. Die Bestimmungs- 
gleichung 


Ag An. 0 


ergibt die Lösungen der ebenen Potentialgleichung für den Kehrwert von f. Der zweite Fall, die 
konstante Schubkraft/Längeneinheit bringt mit 


Re re a EN EB aM SR 5) 
die Differentialgleichung 


und die Lösungen 


RT RE)  nnnl 
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Die verkürzte Differentialgleiehung 
Ergänzt man die Gleichung (7) auf der linken Seite 


2 PP — Syn‘ 1 
@:bu - wD 1E2 BE PESEEuEEEZE (18), 
so tritt an deren Stelle die bedeutend einfachere Form 
F 
a7) = re (19). 
Hieraus folgt 
e AF=T DZ A Eee (20) 
un 
in m F, + F} Pe ee 4 (21). 


F, und ® sind harmonische Funktionen der Ebene und 'F, die Lösung der inhomogenen 
Gleichung (20). 

Für das hyperbolische Keilprofil (10) mit A, = 0 gilt Gl. (19) exakt. Um die verkürzte 
Form auch für andere Profile verwenden zu können, muß man darauf achten, daß der Ergän- 
zungsausdruck (18) vernachlässigbar klein bleibt. Das heißt aber, daß entweder die Komponente 
co, nicht groß sein darf, also 


Fan Ü su 2 A m nee a a a 
gesetzt wird, oder in einem gewissen Profilbereich gilt 
2 u—h “l 

Bu El re ee 


Man ersetzt also das exakte Profil in dem begrenzten Bereich durch die Hyperbel (10), wo- 
bei die Konstanten A, und B, zur Verfügung stehen. 


Lösungen mit dem Produktansatz 


In kartesischen Koordinaten führt das Profil 


er a ee 
mit dem Ansatz za 
F= elttu:y re 


auf die Bestimmungsgleichung für A und 


“+22 + WM—2[K, (HD) HK (Ara +ud)] + -»- 


P4 4 D 1 2 9R 
ne (RR RE) (A) ( + = [Kt —2-K,-Kr-A-u+K2.2]—-0 | (6). 
Für /f = e®Y hat H. D. Coswar!) die Lösungen zusammengestellt. 
Bei beliebigem /(y) erhalten wir, wenn anstelle von e“Y allgemeiner g(y) gesetzt wird 


v1 ‘ m ‘ D} 2. a — R 
9 2.4 +@+! 1 - 


hc cin... 


Zum Hyperbelprofil / = 1/y gehört beispielsweise die Lösung 


RUE 


F=[A,-e&* + B,:e*"]-[G,-cosAy+G,- snAy-+-»- 
+ Ca, ’<sinAy- (2-In2Ay + Ci@ßAy)—iIny2Ay}— cosAy- Si@AY)» +". (08 
+ Ca ssinAy- Siü2Ay) +cosAy-{C (2Ay)—Iny2Ayy)] 2; 


In Polarkoordinaten sind verschiedene Untersuchungen an rotierenden Scheiben bekannt 


!) H.D. Conway, Analysis of Plane Stress in Cartesian : 5 
Arch. XXIX (1960),3. ysıs of Plane Stress in Cartesian Coordinates and with Varying Thickness, Ing. 


ge 
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Für das Profil f= 1/a- r? + b hat E. Anams?) einen Lösun ü 
| 1 ba gsweg für den Fall angegeben, 
daß die Belastungen am äußeren Scheibenrand von r und p abhängen, nämlich in der Form 


Oro) =2&Rr .cosK.nop 
3 PS (207: 
Troy) = 21x -snK.ny 


Schließlich wollen wir von dem allgemeinen Fall, daß auch / von r und o abhängt, al 5 
spiel das Keilprofil betrachten. Mit / Pal 27 


=y—r-sin a ee te ee en 
Re i=y% p (30) 
re). 0. ee EN REEL) 
erhalten wir die an der Spitze eines linearen Keiles angreifende Einzellast (Bild 6). 
Fe Nach Gleichung (7) resultieren aus der Lösung für g(p) die Spannungen: 
| 1 [ 1\ I 
D =, . . ——|e 1 
&y)\ ee A cos (V2 (1 = e) + B-sn(\2- (1) 9). 
en (32). 
= 0 
I 
Ss me n 
\\ > > Immun 
Ne 
=> =. e 
\ —> D ] 
—, / z 
| U, 
o et A 
7 
Ü 
ee 
OEL eu DI 
2 = 
tr 


Bild 6 Bild 7 
Nach der verkürzten Gleichung (19) ergibt sich vergleichsweise 


1 
= 7:[A:c0s2p + B-sin2p], 


RN 


=, 
Pe) 

Wenn P = P, und y, = 0 gesetzt werden, dann befindet sich die maximale Spannung an 
der Stelle 9 = 0 (y = 0), und es gilt für kleine Winkel », die Näherung (22). Die auf osx,max be- 
zogenen Spannungen Oyerk,max und oc max (= (P/F) + (M/W)) wurden über y, aufgetragen (Bild 7). 
Es zeigt sich, daß die Spannung nach der verkürzten Differentialgleichung bis y = 0,5 um nicht 
mehr als 5% vom genauen Wert abweicht. 


Manuskripteingang: 24. 10. 1960 . 
Anschrift: Prof. Dr.-Ing. W. VockE£, Dresden A 16, Bundschuhstr. 1 


2) E. Apams, Der ebene Spannungszustand in einer durch beliebige Randlasten beanspruchten rotieren- 
den Vollscheibe mit dem Stärkeprofil h(r) = 1/(a r? + b), Z. Flugwiss. 5 (1957), H. 11, und Z. Flugwiss. 6 (1958), 
H. 3. 
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KLEINE MITTEILUNGEN 


Bemerkungen zum Iterationsverfahren 
von Schulz zur Matrixinversion 


Faßt man die Inverse einer nichtsingulären, r-rei- 
higen Matrix A als Nullstelle des durch 


DE a EEE? (1) 


definierten Operators 7’ auf und zieht zu ihrer Berech- 
nung das Newronsche Iterationsverfahren [3]}) 


Kar = Xn—{Tany!TXn (n=0,1,2,...) (2) 
heran, so ergibt sich wegen 
TH =X1HX" .....(a) 
I (3b) 


die bekannte Iterationsvorschrift von G. ScHuLz [8], 


[6], [5], [1] 


[Ans = An? B— A X) (n=0, De] (4), 


zu der kürzlich W. Dück [4] eine Fehlerabschätzung 
angegeben hat. In dieser Mitteilung soll gezeigt wer- 
den, daß zu ihrem Beweis (vgl. [4]) auf die i.a. nur 
unter erheblichem Arbeitsaufwand nachprüfbare 
Voraussetzung det X, = 0 verzichtet werden kann; 
ferner, daß man bei der Inversion von ‚‚Matrizen mono- 
toner Art‘‘ unter geeigneten Voraussetzungen auf eine 
monoton wachsende Iterationsfolge geführt wird. Ab- 
schließend werden Ergebnisse der Anwendung des 
„verbesserten NeEwTonschen Verfahrens‘ zusammen- 
gestellt. 


1. Die Fehlerabschätzung 
Satz. Unter der Voraussetzung?) 
#—AX|<1. (5a) 


gilt für die Elemente X, der durch das ScHuLzsche 
Iterationsverfahren (4) erzeugten Folge neben 


im Aue Ar 
[also lim ||E — A X„|| = 0, d.h. lim X, = 4-1 (6)] 
Nn—>XO Nn—>X 


die Fehlerabschätzung 


| : IB—AXn| , 
Ar— Inu = 8a x] A — Zul m. 


Beweis. Die quadratische Konvergenz (5b) folgt 
bekanntlich aus 


E—AX,=(E —A An-ı" = soo (E—A X)” (8). 
Weiter ist 
NEE EN, 


= (A1— X) A X 
er Xn er Anıı 4A Xn 


: 9) 
= (D— (20 - AXNAZN 
= X, (E— A Xy)? 
= (Kun A 
Pe ee. BEE DE bee 1; 
‘) Ziffern in eckigen Klammern verweisen auf das Literatur- 
verzeichnis. 
2): Die Norm ||4|| einer Matrix A hat folgenden Axiomen [5] zu 
genügen (0 bedeutet die Nullmatrix, e eine Konstante): 
1. |4]|>0 für A+O, 
A|=0 für A=0. 
2. |leAll = lel- |lAll. 
8. IIA+B|| = Al] + I12]|, 
l21l s |lAl|- |12]|. 
r r 
Beispiele: |]A|| = Max & |«;,| ||4|| = Max Z& [a;,|° 
s \ = 4|| =Max 2% |a;,| 
ee Eimt ”% 
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Man kann also abschätzen: 
A — Zul {1 — ||# — 4A Xul|} 


<| rn — nl | E A| - - - - - (10) 
und darf wegen ||E — A Xn|| < 1 dividieren: 
1a — 4 Xı|| 


4 — Knril| = TE ar 


nu —Zu]| N). 
2. Inversion von Matrizen monotoner Art 
Definition (J. Schköper [7]. Eine Ordnungs- 


bezieh 
Bi I=B (12a) 


zwischen zwei. Matrizen A= (aj;)) und B= (br) 
bedeutet: 

Die Indizes j und k sind (in gleicher Weise) je in 
zwei elementefremde Klassen aufgeteilt (von denen 
eine auch leer sein darf), und es gilt 


aje <bjr wenn j und k zur selben Klasse gehören 
(12b), 


ak > bjr wenn j und k zu verschiedenen Klassen 
gehören (12e). 


a 0 a 0 a m. wa 


Definition (L. CoLLarz [2]). 
Eine Matrix A heißt ‚‚von monotoner Art“, wenn aus 


AUSAYV (13a) 


ee er 


folgt 
(13b). 


(Matrizen monotoner Art besitzen [2] eine Inverse 
ABO, 


Satz. A sei eine Matrix monotoner Art; dann gilt 
für die Elemente X, der durch das ScHuzzsche Ite- 
rationsverfahren (4) erzeugten Folge unter den Vor- 
aussetzungen 


0OSETAZ, 


Pr u ee 


und "O8 72 


neben 


OSE—AX, 


die Monotonieaussage 


ie, er © 


Beweisdurch vollständige Induktion. (15) und (16) 
seien für einn=p> 0 richtig [Verankerung durch 
(14)], dann rechnet man leicht nach: 


)OSR-AR, 

0 SIE Ze Ar 
Dosen 

B<2E—AX,, 

ApSKpl2E—ANX,) = Xpu; 
ART re Re 

AXpı SAAT, 

FRE PEN 


3. Das „verbesserte NEwronsche Verfahren“ 
Das „verbesserte Newroxsche Verfahren“ [3] 
Tr + 716% = 0 


: Re n= 2 
Ta + een. Ve 


führt hier?) auf die Iterationsvorschrift: 
X = Ant In (E—A Xn) 
Xp = Xn + X (B—AX,) 
u} Tip HK=—XIıHYTı KA YıRrr img 


Iu=0 1, 20a tie), 
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Unter der Voraussetzung ||E — A X,|| < 1 herrscht 
"kubische Konvergenz f 


B—AXRnı = (EB — AR, : =. .(19), 
_ und es gilt die Fehlerabschätzung 


4 IB —ARn| 
At Kalle ar Anal (20. 

Der für Matrizen monotoner Art beim gewöhnlichen 
Newroxschen Verfahren bewiesene Monotoniesatz 
überträgt sich wörtlich auf das verbesserte NEWTON- 
sche Verfahren. 
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High Order Correct Difference Schemes 
for Multi-Dimensional Anisotropic Diffusion 
Equations*) 

In Reference [1], formulae have been deduced for the 
coeffieients in high order correct difference representa- 
tions for the one-dimensional diffusion equation 


ou ou 


nee | al) 
in the form 
U(x, k) = U(2,0) + $ 423 U(x,0) A;ki (2), 
where Ira 
42° u) = 2 TV) ul), | EN 
Aue) =u@ + h)— 2u@) tu@— | 


and the correction factors A; are determined by the 
compatibility conditions for a truncation error of 
order ß 


r B ö 
“—} AI A (4) 
! _.k=0 
with - ; 
tk EN a j 
a yo _.1ı9+%k = r ne — ET 
S. £ za, 1) (, ar 5 (2 s)! E) h2 ( ) 


The first few of the coefficients A; were given in the 
form of the polynomials in 1/R [1] 


Al, 

A=z ( Ber 

Ai! Rt) (6). 
4,= F ( : E vr um) 

A 2 ( nt . ee: T vor m) 


*) This research was supported by the Air Force Office of Scientific 
Research of the Air Research and Development Command. 
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It is readily shown that the three-dimensional an- 
isotropic diffusion equation 


DU. my OU, - 30u, u 
Fra ir ker an) 
has the corresponding difference formula 
Ulz, y, 2, k) 


a a 
= 2 BR Br Ri, RaRzA,_e_ıR1 
SA, IA uam) = (S)% 
where - 
ak ßk yk 
Te za 
and 


Ar eh A;_x_ı(Rı) Ar(R,) 4ı($,) (10) 


and the A; are the coeffieients (6) of the one-dimen- 
sional case. 

In order to prove (10), one only has to compare the 
compatibility equations for (7) 


Bar nee, Be FE Fi 
(.—r —s)!r!s! 


j—k 
j=0 k=0 en Aa duhs 


a a N ee 


(1) 


=0,1,..,n, 1 les, S=E0l,n 


with the compatibility equations (4) for the one-di- 
mensional case. 

The formula (10) follows then immediately from the 
rules for algebraic operations on sequences. 


Beference 
[1] I. R.M. RADok and R. F. MERRIL, Numerical Solution of Boun- 
dary Value Problems. ZAMM 40 (1960). 


Verfasser: R. HERMAN, J.R.M.RaApox, and K. Wang, 
Polytechnie Institute of Brooklyn Com- 
puting Laboratory, 333 Jay Street, Brook- 
lyn 1, New York 


Genäherte Integration in der Nähe eines ein- 
fachen Pols 


Zur genäherten Berechnung des Integrals 
b 
NER N) 
a 


verwendet man meist Summenausdrücke der Form 


N 
J, = 3 f(x) 
k=0 


. (2), 


wobei die Konstanten c7; durch Integration der Stütz- 
polynome L;(x) der LaGranGzschen Interpolations- 
formel 


fa) = 2x) Ka) + Ra) ©. 8) 
gefunden werden: 
b 
Geil rlo)de or kahl: 


a 


Es wird nun angenommen, daß f(x) in x = p einen 
einfachen Pol besitzt. Liegt dieser in der Nähe des 
Integrations- oder des Interpolationsintervalls, so 
führt die genäherte Integration nach (2), (4) natur- 
gemäß zu schlechten Resultaten bzw. zwingt zur Be- 
rechnung einer großen Zahl von Funktionswerten. 

Nun liegt esnahe, die LaGRANGEsche Interpolations- 
formel nicht auf f(x) selbst, sondern auf die Funktion 


gaa)=P— MR)... ... 0) 
anzuwenden, welche den störenden Pol nicht besitzt, 


at) = & In) gen) + Re). . > 


18* 
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Als Näherung für das gesuchte Integral 
b 
I = | ES (7) 
p—x 
a 
erhält man jetzt e 
nr [Ran - - 8) 
ve I 2T* 
a 
oder R 
J= Zrf) een (9) 
5 k=0 
mit s 
P% 
+=| — Io)da. . . « (10) 
a 


Es wird nun gezeigt, daß zwischen den beiden Nä- 
herungen J, und J, eine einfache Beziehung besteht, 
die sich oft mit Vorteil ausnutzen läßt. Es ist 


n 
Ja Jı = 2 —c)far) - - . (A). 
Für die Differenz c$ — c, erhält man aus (4) und (10) 
b 
x ——% 
‘—-%,= j en; Lr(z) de . (12) 


a 
la) = ( —%) (a — a) (am). - (13) 

gilt bekanntlich für die Stützpolynome 
px) 


Mit 


T = — — ı  — — (+ 0%; k=0,1,...,n) (14). 
re TC Baal 
Damit ergibt sich 
= Se RU, 
N | (pP — x) p’(«x) z 
a 
und dies wird in (11) eingesetzt, 
pl®) ”, Mar) 
Je dı = dx» - HERE): 
er 375 Dr) . 


[77 
Die Summe ist aber gerade die Steigung (dividierte 
Differenz) n-ter Ordnung der Funktion f(x), die mit 
den Stützstellen x, bis x, gebildet wird, 


A Earl 
b 


(= Ike dr. 
p—x 


[27 

Will man also bei der genäherten Integration in der 
Nähe eines Pols die Formel (9) an Stelle von (2) ver- 
wenden, so braucht man gar nicht die n -- 1 Kon- 
stanten cf nach (10) zu ermitteln. Man wird vielmehr 
J, nach der üblichen Formel (2) berechnen und daran 
die Korrektur (17a, b) anbringen. Der Faktor C hängt 
vom Ort des Pols, von den Stützstellen und vom Inte- 
grationsintervall, aber nicht von f(x) ab. 


. (17a), 


. (17b). 


I——— 
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Hat man es mit äquidistanten Stützstellen 
en hi (ka 8) 
zu tun, so gilt bekanntlich 
OR 
[%o: % +» +» In] = m 


wobei ö,, die n-te Differenz der (geordneten) Funk- 
tionswerte f(z;) ist, und nach Einführung von 
t= (c — x,)/h und einigen Umformungen erhält man 


für (17a, b) 

ET ED, (208), 

B n 

t dt 
Da (ra (20b), 

& . 

d—%, bb P—u 

=, =, - EZ (200). 


Die Berücksichtigung mehr als eines Pols von f&) 
ist in analoger Weise möglich, jedoch hat die Korrek- 
tur J, — J, dann nicht mehr die einfache Form (17a,b), 
Es läßt sich zeigen, daß (bei Lasrangzscher Inter- 
polation) die Funktion g(x) aus (5) die einzige ist, die 
zu einer Darstellung der Form (17a) führt. 

Beispiel: Sımpsonsche Formel 

Hier gilt die Stützstellenverteilung (18) mit n = 2, 

und esistt«=0,ß8=2. Damit wird 


N d 
t £ 
»=3% (3), 
0 


ing y—2 
u a a le ze u & 
In der beigegebenen Tabelle sind für A = 0,05 und 
verschiedene Integralgrenzen die Näherungen J, und 
J, des Integrals 
b 


J= | Fe — — artanh b — artanh a 
De 


@ 


angegeben. Die drittletzte Spalte gibt die Korrek- 
turen, die von J, zu J, führen. In der letzten Spalte 
steht der Restfehler; die Unregelmäßigkeiten im Gang 
der Werte sind darauf zurückzuführen, daß J, und 
Jg — J, je mit einer Unsicherheit von 1078 behaftet 
sind. Die letzte Zeile soll zeigen, daß auch bei weiterer 
Annäherung an den Pol eine Verkleinerung von h 
meist nicht nötig ist. Der Restfehler J, — J nimmt 
für y— 2, d.h. wenn sich die obere Grenze des Inte- 


grals dem Pol beliebig nähert, nur bis zu 2, _— um 


— — 0,000 005 62 ab. 39 19 

Wie man vorzugehen hat, wenn der Pol im Inte- 
grationsintervall liegt und J als Cavc#vscher Haupt- 
wert gedeutet wird, soll in einer anderen Veröffent- 
lichung gezeigt werden. Der Fall, daß man an Stelle 
von g(2) eine allgemeinere Funktion s(x) - f(x) inter- 
poliert, wird später ausführlicher erörtert. 


Verfasser: Prof. Dr.-Ing. G. Orırz, Dresden A 27, 
Zeunerstr. 88 


ab L A Die | 8 (7, 108| J, (II): 108 
PERS VOR Wenn 2 ir ee ie 
) 0,1 | 0,100 33544 | 20 | —0,01849 0,005 09 ige 9 | 0,100 335 35 ) 
0,1 | 0,2 | 0,102397 32 | 18 | —-0,023 10 0,005 73 . 13 | 0,102 397 19 9 
0,2 | 0,3 | 0,106 78724 | 16 | —-0,029 69 0,007 23 _ 21 0,106 787 03 ee” 
0,3 | 0,4 | 0,11412970 | 14 | —-0,039 55 0,010 18 e 40 | 0,11419390 | — 3 
0,4 | 0,5 | 0,125 658 06 12 —0,055 29 0,015 97 — 88 | 0,125 657 18 — 3 
0,5 | 0,6 | 0,143 843 34 10 | 0,082 74 0,028 45 — 235 | 0,143 840 99 —.5 
0,6 | 0,7 | 0,17416152 | 8 | --0,137 26 0,060 08 _ 825 | 0,17415827 |  — 8 
07.2150,8 | 0,231 356 99 6 0 27188 0,167 13 — 45 35 | 0,231 311 64 | — 12 
0,8 | 0,9 | 0,37425584 | 4 | —0,778 36 0,833 73 — 648 94 | 0,373 60690 | — 30 
0,89 | 0,99 | 1,490 425 37 | 2,2 | — 7,014 4491 | 37,879 1305 | —265 701 23 | 1.224 724 14 — 240 


N 5 
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Beitrag zum schwach gekrümmten Stab unter 
 Außendruck 


Die Tangentialverschiebung v® der Querschnitte 
eines schwach und gleichmäßig gekrümmten Stabes, 
der durch einen konstanten Flüssigkeitsdruck p ent- 
sprechend Bild 1 belastet wird, folgt bei konstanter 


Biegesteifigkeit EJ des Stabes der Differentialglei- 


chung 

dev div dv 

di u 
wenn die Dehnung der Stabachse vernachlässigt und 
mit ausreichender Genauigkeit e = 0 gesetzt werden 
kann, wie dies z. B. beim Knicken des Bogenträgers 
‘ während des Überganges von der stabilen zur in- 


differenten Gleichgewichtsform möglich ist [1]. In 
(1) bedeutet 


meell rk - (1), 


Kal 


sie e Zara. 


Bild 1 


Die vorausgesetzte Dehnungslosigkeit der Stabachse 
schränkt den Anwendungsbereich von (1) ein. Bei 
manchen Problemen läßt sich (1) trotz zunehmenden 
Einflusses der Dehnung noch mit gutem Erfolg an- 
wenden [2]. Bei anderen Problemen wird das Ergeb- 
nis entscheidend von der Dehnung beeinflußt. Sie 
werden daher mit anderen Mitteln gelöst [3], [4]. 


Im folgenden wird eine strenge und geschlossene 
Lösung für den schwach gekrümmten Stab unter 
konstantem Flüssigkeitsdruck nach Bild 1 angegeben. 
Anstelle der Verschiebungsgrößen wird die Krüm- 
mung als Veränderliche benutzt. Dasselbe Ergebnis 
läßt sich jedoch auch mit der Dehnung erzielen. 


Bild 2 


Nach Bild 2 ergibt sich die Dehnung des Stabele- 
mentes zu 


% 


„uU+tr) 0). 
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Es zählt e > 0 bei einer Verkürzung des Elementes. 
Striche kennzeichnen Ableitungen nach g. 


Die Normalkraft im Stab wird mit (3) 
N=BFRe=Eri—"(+y|..(. 
Als Druckkraft ist N >d. 
Das Biegungsmoment ist 
IMS EN (Ze (5). 


Es ist positiv, wenn es die Krümmung des Stabes 
vergrößert. 


N a ‚Eldprdıy) 


Bild 3 


Nach Bild 3 ergibt sich aus dem Gleichgewicht der 
am Stabelement angreifenden Kräfte gegen Drehung 


um O0 


BI ee de m 


Das Gleichgewicht der Kräfte in radialer Richtung 
liefert 


NU+mW-Q=-Zu4m) ...(M 
und in tangentialer Richtung: 


N +QU+W)=0 . (8). 


Mit (4) und (5) wird unter Berücksichtigung von 
” = J/F aus (6) 


mit der Integrationskonstanten cı folgt. 
Geht man mit diesem Ergebnis unter Verwendung 
der Abkürzung 


ie 


in (7) ein und setzt dort weiter Q’ nach (8) und N 
nach (4), so wird zunächst 


h’ AT, p 
Te =nph 


‚2 
und daraus mit W = —. -x x' nach (11) weiter 
0 


Ey a 
Une) ge 


DE yg / 
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Führt man die Differentiation in der vorstehenden 
Gleichung aus und setzt zur Abkürzung darin 


2 
I(x) = EB ’ e 


EJ - 


dann ergibt sich für das Problem die folgende Diffe- 
rentialgleichung: 


| a 0. A: 


Zur Lösung führen wir in (13) die neue Veränder- 
liche z(x) = x’? ein, setzen also 


= (x)? 


„Me 
= — gg 


und 
—1/2 


„ 


n' = a (x) .2(r) we = 5 2(x) , 


und erhalten so die lineare Differentialgleichung 
(x) +2 H)2() H29)=0.. - - (14). 


Die Punkte bedeuten Ableitungen, nach x. Die all- 
gemeine Lösung von (14) lautet 


(x) == Fe) dx fi [es 9 2 [g9(x) & BIP IC) dx e dx] (15). 


Ti2 lo 


In dem Integral [ f(x) dx = j' h(x) 


Zähler gleich der negativen Ableitung des Nenners 
und daher 


- d« ist der 


S I) da = —luh. : 

Somit wird aus 
ESF) de _ Anh _ 7, 
e2Sfddr _ 52h _n2, 


und damit unter Berücksichtigung von (11) und (12) 
aus (15): 


z(x) = h? (« dr — Hr - Le - «) (16). 


Hieraus ergibt sich mit z(x) = .x'? als Lösung von (13) 


(17). 

Die Funktion = g(x) nach (17) und ihre Um- 

kehrung x = x(p) stellen die gesuchte Lösung dar, 

auch wenn sich das Integral (17) nicht durch elemen- 
tare Funktionen in endlicher Zahl darstellen läßt. 
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Zur ebenen Potentialströmung infolge Bewe- 
gung einer geknickten Platte 


Bei der Steuerung von Flugkörpern und Schiffen 
werden Bauteile verwendet, die man — sofern sie 
dünn sind — als geknickte Platten bezeichnen kann. 
Wir wollen nun das durch Ruderbewegungen, z. B, 
Schwingungen, in einer an sich ruhenden Flüssigkeit 


u 


induzierte instationäre Geschwindigkeitsfeld ermit- 
teln. Dazu ziehen wir die Methode der konformen 
Abbildung heran. i 
Bei zeitlicher Änderung der Ruderlage entstehe in 
der ausgangs ruhigen Flüssigkeit eine Strömung mit 
dem komplexen Potential 2, das über die beiden 
folgenden Bedingungen festgelegt wird: 
I) 2 ist im Strömungsgebiet regulär, 
II) Längs der Platte ist die Normalkomponente der 
Konturgeschwindigkeit gleich der Strömungs- 
geschwindigkeit, d. h. 


&2 = 
Be 
dz /n öt /n 
Dabei wird durch 


= |l+r)) se 


die Platte der z-Ebene auf den Einheitskreis der Z- 
Ebene abgebildet. A(t) mißt im Wesentlichen den 
Knickwinkel. Wir bringen nun nach v. Gorvr [1] die 
Randbedingung (1) in eine für uns günstigere Form, 
indem wir sie als Randwerte der analytischen Funk- 
tion (2’/z’)* deuten und von beiden Seiten der Glei- 
chung (1) mittels Normalenvektor auf 
’%* j 2.20% 
Er el 2 0 
L Ü 
übergehen. Diese Funktionalgleichung legt zusam- 
men mit Bedingung I) das 2’ fest zu: 


2' = R(2*z’) + sr 2) = re ha 
Die bisher benutzte Symbolik bezeichnet: 


a vor =Ag); 


+0 —1 +0 
IM) = zu" =- zn" + gm”. 
m Tree ee 
k{) S(f) 
Der Ausdruck (4) ist demnach in £ analytisch. 


Das Geschwindigkeitsfeld infolge des Ruderaus- 
schlags ist gegeben durch (2), (4) und 


Um eine Steuerwirkung zu beobachten, muß man 
dieses % nun — unter Berücksichtigung der Kurtra- 
schen Abflußbedingung — der Umströmung des Ru- 
ders linear überlagern. 


Ast 
I 
Iı\ 
I\ 
| As3 
IIı\ 
Jam 
14 
..r 


Bild 1, Sonderfälle geknickter Platten 


® 


Die ausgezogenen Linien von Bild 1 zeigen die 
Ruderform z. Z. t = 0, die gestrichelten eine dazu be- 
nachbarte zur Zeit t + 0. 


Aus der Abbildungsfunktion (2) ergibt sich: 


LAN 1 TAsmE 
lat) ne +2): 


& 
(6). 


ö 
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 Buchbesprechungen 


Die Aufspaltung in R und 8 wird füri=1 ;3/2;2 eine große praktische Bedeutung zukommt. Wir er- 
besonders einfach, wenn man beachtet, daß rechnen aus (6) das Geschwindigkeitsfeld: 


Unsere Betrachtungen seien jedoch auf die oben- 
genannten Sonderfälle beschränkt, da gerade ihnen 


> 


Bild 2. Elliptisches und JOoUROWSKI-Ruder 


1 1 EA a Ba me el bir 

n(e+2)-ma+ Ss marg+mlııt) A lazı E e+7 a = 

ah = 0 Zpae em, rar 
| Bei allen andern A+3/2,1<A<2 bereitet die Ale el 7) "Tre 
E Zerlegung einige Schwierigkeiten. Schreiben wir i et 
u ee Re an, 
# e+sy-adren, en 
| ee+ 7) je 
8 dann ist auch hier unter Hinzunahme von 
3 +00 ö 1 ( 112, (+02 
; sin & ız n+1 8% — — |—[(C + z) In 
ed sec Tanlektrl Te 
4 —o h 10 1 
| Dane} a (ek up an 
> die gewünschte Aufspaltung von 2*z’ möglich; z. B. = r 3 & 
5 führt diese Darstellung für « = 1/2 auf: 2 | ie 1 
; a E 
2 Fra = 1 = 
B Y= rg [oe tg VE + arc tg E > 2 läßt sich für ein symmetrisches Polynomprofil 
4 A=1 


(Bild 2) genau so einfach hinschreiben, wie man aus 
der zugehörigen Abbildungsfunktion 


N 
BEHZ RTN 


erkennt. Die Druckverteilung am Profilrand findet 
man über die geeignet modifizierte BERNOULLI-Glei- 
chung, [2]. 
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BUCHBESPRECHUNGEN 


E.L. Lehmann, Testing Statistical Hypothe- 
ses. XIII + 3698. New York 1959. John Wiley & 
Sons, Inc. Preis geb. $.11.00. 


Das vorliegende Buch ist aus Vorlesungen des Verf. 
hervorgegangen, die schon früher einem kleineren 
Leserkreis in mimeografierter Form zugänglich waren 
und deren Einfluß in vielen Veröffentlichungen zu 
erkennen war. Ich möchte daher annehmen, daß 
dieses Buch den Charakter eines Standardwerkes für 
die moderne Testtheorie annehmen wird. Das Werk 
zeichnet sich vor vielen anderen dadurch aus, daß 
es in mathematischer Hinsicht keine Kompromisse 
macht. Trotzdem ist es, eine gewisse Höhe des ab- 
strakten Denkens vorausgesetzt, fast voraussetzungs- 
los lesbar. Ein ganzes Kapitel ist den wahrschein- 
lichkeitstheoretischen Grundlagen gewidmet, und in 
einem Anhang ist ein weiterer Teil des benötigten 
mathematischen Werkzeuges dargestellt. Das erste 
Kapitel ist einer Beschreibung der allgemeinen Ent- 
scheidungstheorie gewidmet, und. auf diesem Hinter- 
grund werden dann die Neyman-Pearsonsche 
Testtheorie und ihr Ausbau sowie die duale Theorie 
der Konfidenzbereiche dargestellt. Eingehende Wür- 
digung erfahren die Theorie der trennscharfen und der 
unverfälschten Teste. Vom besonderen Interesse 
scheinen mir die Kapitel über die Invarianztheorie 
und das Minimaxprinzip zu sein. Das letztgenannte 
Kapitel enthält auch eine Darstellung des Theorems 
von Hunt und Stein. Es sind stets viele und wich- 
tige Beispiele in den Text verflochten, welche das 


Buch auch für jenen Leser ansprechend machen, der 
vor allen Dingen an den Anwendungen interessiert 
ist, dabei aber den Dingen auf den Grund zu gehen 
wünscht und nicht nur eine Formelsammlung sucht. 
Das Werk enthält auch für den Kenner interessantes 
und neues Material. Man legt es mit dem Bedauern 
aus der Hand, daß der Autor, der mit diesem Werk 
seine Darstellungsgabe dokumentiert hat, sich nicht 
entschließen konnte, auch noch die Theorie der 
Schätzungen mitaufzunehmen. 


Wien L. SCHMETTERER 


G. Dahlquist, Stability and Error Bounds in 
the Numerical Integration of Ordinary Dif- 
ferential Equations. 85 S. Stockholm 1959. K. 
Tekniska Högskolans Handlingar. 


Das vorliegende Heft enthält Ergebnisse des Ver- 
fassers über die Stabilität von Differenzengleichungen 
zur Lösung von Anfangswertaufgaben bei gewöhn- 
lichen Differentialgleichungen. Anknüpfend an eine 
frühere Arbeit (Math. Scand. 4 (1956), S. 33—53) 
werden Gleichungen der Form 

k 
a (m Mn) = 
zur numerischen Lösung von yfr) = f(x, y), sowie zu 
y'' = f(x, y, y') gehörende Differenzengleichungen be- 


- | 
"trachtet. Den Untersuchungen werden die Begriffe 
instabilen Operators (Nullstellen des Polynoms e(£) 
eo 5 % & sind dem Betrage nach <1,> lund 1, 


ro im letzten Fall das Gleichheitszeichen minde- 
stens einmal angenommen werden muß; dabei werden 
schwach instabile Operatoren also noch stabil ge- 
nannt) zugrunde gelegt. Es zeigt sich, ‚daß die Ap- 
 proximationsordnung p im Vergleich nicht zu groß 
sein darf, genauer: ein Differenzenoperator ist sicher 


stark instabil, wenn p>2 (4 | + E (r + )) 


‚gilt. Hiernach kann z. B. die Sımpsonregel p=4 
k=2,r=1) stabil sein. Wie nachgewiesen wird, 
ist sie insbesondere schwach instabil (8- auch QUADE, 
ZAMM 39 (1959), S. 117—134; dort wird der Begriff 
der Stabilität enger gefaßt, so daß die Sımpsonregel 
instabil ist). 

Neben diesen Abschnitten über starke und schwache 
Instabilität enthält die Arbeit Ausführungen über den 
Verfahrensfehler und über Fehlerschranken bei sta- 
bilen Operatoren sowie eine Zusammenstellung von 
numerischen Beispielen. 

Zusammen mit der genannten Arbeit des Verfassers 
ist dieses Heft für jeden sehr empfehlenswert, der 
sich mit dem behandelten Themenkreis zu befassen 
hat. 


Dresden J. W. SCHMIDT 


P.I. Richards, Manual of Mathematical Phy- 
sies. X + 4868. London 1959. Pergamon Press. 
Preis geb. £ 5.10 =. 


Eine Darstellung, die es sich, wie es im Vorwort 
heißt, zum Ziel gesetzt hat, ‚‚die theoretische Physik 
(so weit sie gesichert ist), ihre Anwendungen und ihr 
mathematisches Rüstzeug in einem einzigen Hand- 
buch von ‚‚vernünftigem‘‘ Umfang (450 8.) zusam- 
menzufassen, ohne an dem logischen Zusammenhang 
und den Grundlagen etwas zu opfern“, muß sich 
naturgemäß auf das Allerwesentlichste beschränken, 
auch dann, wenn der zur Verfügung stehende mathe- 
matische Apparat eine sehr kurze und raumsparende 
Ausdrucksweise gestattet. Auf 250 Seiten wird in 
sechs Kapiteln eine Gesamtübersicht über die Theore- 
tische Physik gegeben. Der Mechanik (einschl. der 
Mechanik der Kontinua) sind 28 Seiten gewidmet, der 
Thermodynamik 30, der elektromagnetischen Theorie 
(einschl. Optik und Informationstheorie) 50, der Re- 
lativität 15, der Quantenmechanik 63 und der Stati- 
stischen Physik 64. Der mathematische Teil, auf dem 
das Schwergewicht dieser Besprechung liegen möge, 
umfaßt etwa 200 Seiten. Es ist kein Zweifel, daß es 
dem Verfasser durch Verzicht auf Dinge, die in der 
überwiegenden Mehrzahl der schon vorhandenen 
Hand- und Tabellenbücher gefunden werden können, 
gelungen ist, auch die modernen Hilfsmittel der Ma- 
thematik zu berücksichtigen. Jedoch bleibt die Dar- 
stellung auch hier äußerst knapp. Sie ist in 16 Kapitel 
eingeteilt. Aus der Algebra (12 Seiten) werden einige 
wenige Tatsachen über Polynome und rationale Funk- 
tionen und ihre Anwendung bei Interpolationen, über 
Wurzelabschätzungen und einige endliche Summen 
mit Binomialkoeffizienten angeführt. Die Hinweise 
auf die Nomographie dürften in dieser Kürze wenig 
Zweck haben. Die Differential- und Integralrech- 
nung (17 8.) dringt bis zu den allgemeineren Integral- 
begriffen (Lebesgue, Stieltjes) und zur Theorie 
der Distributionen vor. Außerdem werden hier wie 
im folgenden Kapitel über unendliche Reihen (8 8.) 
asymptotische Entwicklungen besonders berück- 
sichtigt. Esfolgen Kapitel überVektoranalysis (11 $.), 
über Determinanten und Matrizen (14 $., mit einem 
kurzen Hinweis auf den Graphenbegriff), über die 
Funktionentheorie (7(!)S.) und über Integraltrans- 


eines stabilen, eines stark instabilen und eines schwach 


und partiell 
etwas 


TO l 
ee auf den n- onalen 
und im dee ee damit sr! 
und Operatorbegriff eingegangen (14 8.). 
wertproblem und der Störungstheorie sind 
Wahrscheinlichkeits- und der Spieltheorie 14, deı 
allgemeinen Tensoranalysis 11 und schließlich der 
Gruppentheorie (einschl. Darstellungsbegriff) 13 Sei- 
ten gewidmet. An einigen Stellen werden Fragen der 
numerischen Behandlung der anfallenden Aufgaben 
gestreift, aber auch wirklich nur gestreift. Da eine 
Vollständigkeit in einem Rahmen, wie ihn sich der 
Verfasser gesteckt hat, nicht erreichbar ist, wäre es 
unbedingt erforderlich gewesen, den Leser auf ein- 
schlägige Lehrbücher und andere ausführlichere Be- 
handlungen der angeschnittenen Probleme hinzuwei- 
sen. Leider geschieht dies nur äußerst selten und in 
völlig unzureichendem Maße. Auch Namen sollten 
in noch größerem Umfang genannt werden. (Ein 
Beispiel für viele: im Kapitel über die numerische 
Lösung von Variationsproblemen wird ein Ansatz 
gemacht, der überall als Ritz-Ansatz bekannt ist, 
aber der Name Ritz wird nicht erwähnt.) An sich ist 
dieses Handbuch gut lesbar geschrieben, jedoch würde 
es zweifellos durch Hinzufügung solcher Angaben über 
die Literatur beträchtlich an Wert gewinnen. 


Dresden H. HEINRICH 


R.Sauer, Ingenieur-Mathematik. Erster Band. 
VIII + 304 S. m. 178 Abb. Berlin/Göttingen/Heidel- 
berg 1959. Springer-Verlag. Preis geb. 24,— DM. 


Man muß es dankbar begrüßen, wenn sich ein 
Hochschullehrer, der in so hohem Maße wie der Ver- 
fasser im Mittelpunkt des wissenschaftlichen Lebens 
steht, der Aufgabe unterzieht, ein mathematisches 
Lehrbuch für Ingenieurstudenten zu schreiben. Nicht 
nur an der textlichen Formulierung, sondern vor allen 
auch aus dem didaktischen Aufbau spricht eine reiche 
pädagogische Erfahrung. Jede Frage wird erst dort 
aufgeworfen, wo sie der Lernende erfassen kann, und | 
ihr wird gerade so weit nachgegangen, wie es für die | 
Anwendungen oder für das Verständnis des Folgen- | 
den notwendig ist. Der Stoff ist umfassender als aus 
den Überschriften der beiden Kapitelhervorgeht Das 
erste Kapitel ‚‚Differential- und Integralrechnung für 
Funktionen von einer Veränderlichen‘‘ enthält auch 
alles Notwendige aus der analytischen Geometrie der 
Ebene, insbesondere der Kegelschnitte. Es bringt 
neben dem Grenzwertbegriff sofort auch die Grund- 
tatsachen der Reihenlehre, um diese vor allem auf 
die Taylorreihen anwenden zu können. Im zweiten 
Kapitel „‚Differential- und Integralrechnung für Funk- 
tionen von mehreren Veränderlichen‘‘ werden außer 
den im engeren Sinne hierzu gehörigen Sätzen und 
Begriffen auch die Vektorrechnung, lineare Glei- 
chungssysteme (ohne Matrizen) und Determinanten 
und die analytische Geometrie des Raumes, insbe- 
sondere die Flächen 2. Ordnung behandelt. Mit Rück- 
sicht auf die Anwendungen, beispielsweise auf die 
Getriebelehre, wird in beiden Kapiteln großer: Wert 
auf die Differentialgeometrie von Kurven und Flä- 
chen gelegt. In einem Maße, wie es bisher in mathe- 
matischen Lehrbüchern für Ingenieure wohl nicht der 
Fall war, wird auf numerische und graphische Me- 
thoden und damit zusammenhängende Fragen, z.B. die 
Genauigkeitsfrage bei Näherungen, eingegangen. Der 
Verfasser verzichtet nicht auf mathematische Strenge 


(eine Reihe von Beweisen, die den Fluß der Dar- 
stellung hemmen würden, findet im Anhang seinen 
. Platz), aber er verzichtet auf vielen unnötigen Ballast 
und ist nur dadurch und durch die Klarheit seiner 
Ausführungen in der Lage, den Ingenieurstudenten 
ein mathematisches Lehrbuch in die Hand zu geben, 
das einen nicht zu großen Umfang hat und doch das 
Wesentliche enthält. Dieses Lehrbuch wird sich sicher 
bald einer großen Verbreitung erfreuen. 


Dresden H. HEInRıcH 


J. Horn und H, Wittich, Gewöhnliche Diffe- 
rentialgleichungen. Sechste Aufl. 275 S.m. 10 Abb. 
Berlin 1960. Walter de Gruyter & Co. Preis DM 32,— 


Als ausgezeichnet lesbare und leichtverständliche 
Einführung in die Theorie der gewöhnlichen Diffe- 
rentialgleichungen hat sich das Buch von J. HoRN 
schon immer großer Beliebtheit erfreut. Durch die 
neue Bearbeitung, die es jetzt durch H. WırTıca 
erfahren hat, hat sich der Charakter des Buches nicht 
geändert. Es baut auf relativ bescheidenen Vorkennt- 
nissen der Analysis auf und führt von den elemen- 
taren Integrationsmethoden aus allmählich in tiefer 
liegende Fragen der Theorie ein. Dabei wird auch 
auf die numerische und graphische Behandlung, wenn 
auch nur in den ersten Anfängen, und auf Anwen- 
dungen auf schwingende Systeme eingegangen. Da 
sich die Darstellung der speziellen Funktionen, die 
gewöhnlichen Differentialgleichungen genügen, erst 
durch eine Betrachtung im Komplexen richtig ab- 
runden läßt, ist der Theorie der gewöhnlichen Diffe- 
rentialgleichungen im Komplexen wesentlich mehr 
Platz als in den früheren Auflagen eingeräumt. In 
dem Kapitel über spezielle lineare Differentialglei- 
chungen zweiter Ordnung werden Integraldarstel- 
lungen und asymptotische Entwicklungen der Lö- 
sungen behandelt. Wertvoll, auch vom Standpunkt der 
Anwendungen, sind die Ausführungen über Diffe- 
rentialgleichungen mit periodischen Koeffizienten und 
über die Abhängigkeit der Lösungen von Parametern 
und Anfangswerten. Ein Kapitel führt in einige 
Grundbegriffe und Grundtatsachen der qualitativen 
Theorie ein, die, wie z. B. die asymptotische Stabilität, 
für die mathematische Behandlung von Regelungs- 
vorgängen von großer Bedeutung sind. 


Das Buch kann nicht nur Studierenden der Mathe- 
matik als Einführung, sondern auch Physikern und 
Ingenieuren, die in erster Linie an den Anwendungen 
interessiert sind, als Lehrbuch empfohlen werden. 


Dresden H. HEINRICH 


W. R. Hawthorne, W. T. Olson, Design and 
Performance of Gas Turbine Power Plants. 
XIII + 5638. Princeton, N.J., 1960. Princeton 
University Press. Preis geb. $ 15.00. 


Innerhalb der bekannten, von Princeton Univer- 
sity Press herausgegebenen zwölfbändigen Schriften- 
reihe High Speed Aerodynamics and Jet Propulsion, 
welche die für die Luftfahrt notwendigen Grund- 
wissenschaften und ihre Anwendung auf Zelle und 
Triebwerk zum Inhalt hat, nehmen die Bände X 
und XI insofern eine Sonderstellung ein, als sie die 
wichtigsten Baugruppen eines Gasturbinen-Trieb- 
werkes und die in ihnen ablaufenden Vorgänge ein- 
gehend behandeln und damit besonders eng mitein- 
ander verbunden sind. Während jedoch Band X den 
thermodynamischen und strömungstechnischen Fra- 
gen von Kompressor und Turbine gewidmet ist, ge- 
währt der vorliegende, von mehreren Autoren mit- 
verfaßte XI. Band, für dessen Redaktion im ganzen 
W. R. HAWTHoRNE und W. T. Orsox verantwortlich 
zeichnen, einen Einblick in das komplexe Gebiet der 
Brennstoffverbrennung. Auch werden das Verhalten 


Buchbesprechungen 


der Werkstoffe bei hohen Temperaturen und die 
Schaufelschwingungen sowie die Thermodynamik und 
das Betriebsverhalten von Gasturbinen-Triebwerken 
einer Betrachtung unterzogen. 


Der erste Teil enthält einen gedrängten Überblick 
über die Probleme, die in den folgenden Teilen be- 
handelt werden. Am umfangreichsten ist der zweite 
Teil (über 60% des Buchumfanges). Er ist den Fra- 
gen der Brennkammer-Entwicklung gewidmet und 
schließt mit einem Abschnitt, in dem neben der ge- 
schichtlichen Betrachtung über Brennkammerent- 
wicklung das Betriebsverhalten dieses Aggregates 
und ausführliche Hinweise für den Entwurf neuer 
Brennkammern besonders hervorzuheben sind. 


Der dritte Teil (etwa 20% des Buchumfanges) ist 
den Werkstofffragen und den durch aerodynamische 
Kräfte angefachten Schaufelschwingungen, der vierte 
und letzte Teil der Thermodynamik der Gasturbinen- 
Triebwerke und ihrem Betriebsverhalten gewidmet. 
Dabei werden die thermodynamischen Vorgänge in 
Strahlturbinen, Propellerturbinen und Zweistrom- 
turbinen getrennt betrachtet. 


Wenn auch die Vielfältigkeit des behandelten 
Stoffes trotz des gewählten Buchumfanges eine er- 
schöpfende Behandlung einzelner Fragen stellen- 
weise nicht zuließ, so hat dadurch die Darstellungs- 
art an Straffheit gewonnen, und der interessierte 
Leser wird trotzdem durch sorgfältig zusammen- 
gestellte Literatur am Ende jedes größeren Abschnit- 
tes mehr als entschädigt. 


Im ganzen gesehen stellt das vorliegende Buch 
eine wertvolle Bereicherung der luftfahrttechnischen 
Literatur dar und dürfte von Fachkreisen und 
Studierenden, die mit den Grundwissenschaften der 
Luftfahrt vertraut sind, mit großem Interesse auf- 
genommen werden. Daß das Außere des Bandes den 
Gepflogenheiten des Verlages gemäß ansprechend 
ausfiel, braucht nieht besonders betont zu werden. 


Dresden A. PAWLOWITSCH 


Les Math&matiques de I’Ingenieur. Comptes 
Rendus du Congres International, Mons et Bruxelles 
9.—14. Juin 1958. Herausg. L. Derwidue u. N. 
Forbat. Preis brosch. 300 F. 


Die Veranstalter des Kongresses ‚‚Les Mathemati- 
ques de l’Ingenieur“, der im Juni 1958 in Mons 
(Belgien) stattgefunden hat und in Brüssel abge- 
schlossen worden ist, haben die dort gehaltenen Vor- 
träge herausgegeben und damit der breiten Offent- 
lichkeit zugänglich gemacht. Das Schwergewicht 
liegt auf den Hauptvorträgen, welche die Herren 
Bouzırar (Paris), CoLLATz (Hamburg), DE PossEL 
(Paris), Ecervary (Budapest), FicHERA (Rom), 
FiıscHkr (Berlin), JANssens (Brüssel), KIVELIOVITCH 
(Paris), MAssonner (Lüttich-Liege), Rocarn (Paris), 
SAUER (München) und VERNOTTE (Paris) auf Ein- 
ladung der veranstaltenden Fakultät gehalten haben. 
Gemäß der Thematik des Kongresses beschäftigen 
sie sich mit Problemen der Mathematik, die für die 
wissenschaftlichen Ingenieurarbeit von besonderer 
Bedeutung sind. Dasselbe gilt von den 40 anderen 
Vorträgen, in denen in getrennten Sektionen Algebra, 
Analysis, Nomographie und Statistik, Maschinelles 
Rechnen und verschiedene Anwendungen über Einzel- 
probleme berichtet wurde. Die zu den Vorträgen 
gehörigen Abbildungen sind, soweit sie nicht im 
Text untergebracht werden konnten, in einer bei- 
gefügten Mappe enthalten. Den Abschluß bilden die 
auf der Schlußsitzungin Brüssel gelieferten interessan- 
ten Beiträge zu Fragen der Ingenieurausbildung, ins- 
besondere der mathematischen. Das Erscheinen 
dieses Bandes wird zweifellos nicht nur von den 
Teilnehmern des Kongresses, sondern von einem sehr 
viel größeren Interessentenkreis begrüßt werden. 


Dresden H. HEINRIcH 
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A.0. Gelfond, The Solution of Equatione in 


 Integers. 72 8. Groningen 1960. P. Noordhoff Ltd. 


Preis brosch. f 3,75. 


Die — ganzzahlige — Lösung algebraischer Glei- 
chungen beliebigen Grades in mehr als einer Unbe- 
kannten mit ganzzahligen Koeffizienten gehört zu 
einem der ältesten und schwierigsten Probleme der 
Zahlentheorie. Die ersten Untersuchungen in dieser 
Richtung finden sich bereits im Altertum. Die Neu- 
zeit brachte dann eine Reihe sehr schöner Einzelergeb- 
nisse, die mit den Namen zahlreicher hervorragender 
Mathematiker verknüpft sind. Dennoch ist es bisher 
nicht gelungen, allgemeine Methoden zur Lösung der- 
artiger Gleichungen zu finden, wenn der Grad der 
Gleichung größer als 2 und die Anzahl der Unbekann- 
ten größer als 1 ist. Man denke nur an die berühmt 
gewordene FErMmATsche Vermutung. 


Das vorliegende Büchlein des bekannten sowjeti- 
schen Zahlentheoretikers hat es sich zur Aufgabe ge- 
macht, einen größeren Kreis mathematisch interessier- 
ter Studenten, Lehrer und Ingenieure mit einem Teil 
der bisherigen Ergebnisse und den bei den Beweisen 
verwendeten mathematischen Methoden bekannt zu 
machen. Der Stoff ist so klar und elegant dargestellt, 
daß man das Buch mit großer Freude liest. 


Im einzelnen werden in den 7 Paragraphen des 
Buches behandelt die Gleichungen 1. Grades in einer 
und zwei Unbekannten, Beispiele von Gleichungen 
2. Grades in drei Unbekamten (+ yY= 2%, 
x? +2 y? = 2), Gleichungen 2. Grades in zwei Unbe- 
kannten und schließlich einige Ergebnisse über Glei- 
chungen höheren als 2. Grades in zwei und drei Un- 
bekannten. 


Dresden M. HasseE 


L. Collatz, Differentialgleichungen für In- 
genieure. 2. Aufl. 1978. m. 115 Abb. Stuttgart 
1960, B. G. Teubner Verlagsgesellschaft. Preis geb. 
DM 21,60. 


Bei dieser für Ingenieure und für Studierende der 
technischen Wissenschaften bestimmten Einführung 
in das Gebiet der Differentialgleichungen und ihrer 
Anwendungen hat es sich der bekannte Verfasser 
zum Ziel gesetzt, ohne Verletzung der notwendigen 
mathematischen Strenge die anschauliche Seite zu 
betonen und bei den Begriffsbildungen und bei der 
Darstellung der mathematischen Zusammenhänge 
immer wieder Deutungen aus der Mechanik und aus 
anderen Anwendungsgebieten heranzuziehen und auf 
diese Weise die Scheu vor dem Arbeiten mit Differen- 
tialgleichungen überwinden zu helfen. Man darf 
wohl sagen, daß es einen besseren Weg zur Erreichung 
dieses Zieles, als ihn der Verfasser beschritten hat, 
kaum geben dürfte. Die Stoffauswahl, auf die es 
hierbei entscheidend ankommt, ist vom Anfang bis 
zum Ende den Anwendungen angepaßt, und trotz 
des im ganzen geringen Umfangs des Bändchens 
werden viele einfache Probleme aus Physik und 
Technik, die auf Differentialgleichungen führen, als 
Beispiele angeführt und gelöst. Neben den ge- 
schlossenen integrierbaren Typen von gewöhnlichen 
Differentialgleichungen erster Ordnung haben vor 
allem dielinearen Differentialgleichungen zweiter und 
höherer Ordnung eine betonte Berücksichtigung ge- 
funden; vor allem wird auf die dabei auftretenden 
Rand- und insbesondere Eigenwertprobleme großes 
Gewicht gelegt. Damit ist endlich in einer Weise, wie 
es sonst beieinführenden Darstellungen für Ingenieure 
kaum der Fall gewesen ist, diesem so wichtigen Teil- 
gebiet der ihm gebührende Platz eingeräumt. Von 
den speziellen Differentialgleichungen werden die- 


und den hypergeom 
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auch an en erh ; Pr 
eingegangen. Der Anhang ist einigen Nä 
Bean Kam (RungE-KUTTA, Differenzen- 
stellenverfahren, Rırz) gewidmet. ae, ie‘, 
Man darf sicher sein, daß sich diese 
schnell großer Beliebtheit erfreuen wird. 
stattung dieses Bandes I einer neuen Reihe 
Verlages ist ausgezeichnet. u 
H.Hemson 
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Jahnke-Emde-Lösceh, Tafeln Höherer Funk- 
tionen. 6. Aufl. XII + 318 S. m. 189 Abb. Stutt- 
gart 1960. B.G. Teubner Verlagsgesellschaft. Preis 
geb. DM 58,80. j 


Die Tafeln höherer Funktionen von JAHNKE-EMmDE 
haben sich so bewährt und sind in der Welt so bekannt, 
daß es sicher einem echten Bedürfnis entspricht, daß 
sie jetzt neu herausgebracht worden sind. Die vor- 
liegende Neubearbeitung von F. LöscH zeichnet sich 
nicht nur durch große Sorgfalt, Gründlichkeit und 
Zuverlässigkeit aus, sondern sie zeugt auch von einer 
reichen Erfahrung, wie Funktionentafeln zu gestalten 
sind, damit sie dem Benutzer alle erforderlichen Be- 
quemlichkeiten und Erleichterungen bieten. Gerade 
in dieser Hinsicht ist sehr viel getan worden, und fast 
jede Seite verrät ein großes Maß mühevoller Klein- 
arbeit. Alle Tafeln mit doppeltem Eingang für ein 
Argumentsind durch Spaltentafeln mit einem Eingang 
ersetzt worden. Nicht zuletzt dadurch ist es möglich 
geworden, den Funktionswerten die ersten Differenzen 
für die lineare bzw. die zweiten Differenzen für die 
quadratische Interpolation (nach LAPLACE-EVERETT) 
beizufügen, eine sehr wesentliche Verbesserung! Je- 
dem Abschnitt ist über die Definitionen und Bezeich- 
nungen der dort behandelten Funktionen eine text- 
liche Einführung vorangestellt, die in knapper und 
klarer Form das Wesentlichste zusammenfaßt. Der 
dafür beanspruchte Platz ist durch eine etwas ge- 
drängte Gestaltung der Tafeln (ohne Verlust an Über- 
sichtlichkeit!) und teilweise durch Wahl eines größeren 
Argumentschrittes (wo dies ohne Verzicht auf be- 
queme Interpolierbarkeit möglich war) gewonnen 
worden. Die Bezeichnungen der Funktionen sind, wo 
es notwendig war, dem im heutigen Schrifttum vor- 
herrschenden Brauch angepaßt worden, z. B. ber x 
und beix statt Re Ilx yi) und — Im BAR: Vi). Die 
Funktionskurven und vor allem die Reliefdarstellun- 
gen sind wohl im wesentlichen dem alten JAHNKE- 
EmpE entnommen. Auch sonst ist die Stoffauswahl 
bis auf einige durch die heutigen Bedürfnisse bedingte 
Ergänzungen und Kürzungen im großen und ganzen 
dieselbe geblieben, jedoch sind Umstellungen und 
Umgestaltungen vorgenommen worden. Seit der 
letzten Auflage erschienene Arbeiten wurden selbstver- 
verständlich berücksichtigt. Teilweise haben sie 
ganze Abschnitte in Form und Inhalt wesentlich 
beeinflußt. Das reichhaltige Literaturverzeichnis gibt 
auch davon einen Eindruck. 


Es ist kein Zweifel, daß das bewährte Tafelwerk 
durch die Neubearbeitung noch erheblich gewonnen 
hat. Da es zudem vorzüglich ausgestattet ist, wird es 
auch in dieser neuen Gestalt (trotz des Preises!) die 
ihm gebührende weite Verbreitung finden und bald 
ein unentbehrliches Hilfsmittel der angewandten 
Analysis sein. 


Dresden H. HEINRICH 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 
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EINGEGANGENE BÜCHER 
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Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


E. Sehörner, Raumbild-Lehrbuch der dar- 
stellenden Geometrie für Ingenieurschulen. 
153 S. m. 60 Anaglyphen und 228 Zeichnungen. Mün- 
chen 1960. Oldenbourg Verlag. Preis geb. DM 16,—. 


H.-J. Kawalsky, Topologische Räume. 2718. 
m. 2 Abb. Basel 1961. Birkhäuser Verlag. Preis 
geb. sFr. 40.—. 


A. Korganoff, Methodes de Calcul Numeri- 
que. Tome 1: Algebre non lineaire. XXVII + 375 8. 
m. 20 Abb. Paris 1961. Dunod Editeur. Preis geb.. 
58 NF. 


W. Meyer zur Capellen / K.-A. Rischen, Lagen - 
zuordnungen an ebenen Viergelenkgetrieben 
in analytischer Darstellung. Eine Maßsynthese. 
(Forschungsberichte des Landes Nordrhein-Westfalen, 
Nr. 923.) 848. m. 29 Abb. Köln und Opladen 1961. 
Westdeutscher Verlag. Preis brosch. DM 23,20. 


W. Rinow, Die innere Geometrie der metri- 
schen Räume. XVI- 5208. Berlin/Göttingen/ 
Heidelberg 1961. Springer-Verlag. Preis geb. 
DM 83,—. - 


B. Eck, Technische Strömungslehre. Sechste 
'neubearbeitete Aufl. XII +4538. m. 454 Abb. 
Berlin/Göttingen/Heidelberg. Springer-Verlag. Preis 
geb. DM 31,50. 


J. P. Schouten, Operatorenrechnung. Mit An- 
wendungen auf technische Probleme. VIII + 2248. 
m. 128 Abb. Berlin/Göttingen/Heidelberg 1961. 
Springer-Verlag. Preis geb. DM 31,50. 


H. Wolf, Spannungsoptik. Ein Lehr- und Nach- 
schlagebuch für Forschung, Technik und Unterricht. 
XIX + 5828. m. 311 Abb., 1 Tafel u. 3 Tabellen. 
Berlin/Göttingen/Heidelberg 1961. Springer-Verlag. 
Preis geb. DM 66,—. 


S. Krug — P. Stein, Einflußfelder orthogonal 
anisotroper Platten. Influence surfaces of ortho- 
gonal anisotropie plates. In deutscher und englischer 
Sprache. Englische Übersetzung durch H. Juhl. 
VIII + 328. m. 31 Abb. u. 193 Tafeln. Berlin/Göt- 
tingen/Heidelberg 1961. Springer-Verlag. Preis geb. 
DM 67,50. 


L. Prandtl, Gesammelte Abhandlungen zur 
angewandten Mechanik, Hydro- und Aero- 
dynamik. In 3Teilen mit einem Porträt und 
1003 Abb. Teill: XX+S. 1—574, Teil2: XIV + 
S. 575—1070, Teil 3: XIV + 8. 1071—1620. Berlin/ 
Göttingen/Heidelberg 1961. Springer-Verlag. Preis 
brosch. Teil 1-3 DM 273,—. 


R. Sauer, Ingenieur-Mathematik, Band 2: 
Differentialgleichungen und Funktionentheorie. XII 
+ 1808. m. 95 Abb. Berlin/Göttingen/Heidelberg 
1961. Springer-Verlag. Preis geb. DM 18,—. 


F. Hohenberg, Konstruktive Geometrie in der 
Technik. Zweite Aufl. IX + 3198. m. 459 Abb. 
Wien 1961. Springer-Verlag. Preis geb. DM 29,—. 


Handbuch der Physik, Band XI/Teil 1: 
Akustik I. VI+ 4438. m. 183 Abb. Berlin/Göttin- 
gen/Heidelberg 1961. Springer-Verlag. Preis geb. 
DM 128,—. 


M. Neumann, Zur quantitativen Auswertung 
von Schlierenbildernrotationssymmetrischer 
Strömungen. (Forschungsberichte des Landes 
Nordrhein-Westfalen, Nr. 785.) 528. m. 6Tab. u. 
6 Abb. Köln und Opladen 1961. Westdeutscher Ver- 
lag. Preis brosch. DM 15,30. 


Proceedings of the Fifth Congress on Theo- 
retical and Applied Mechanics and the Sym- 
posium on Non-linear Physical Problems, 
Roorkee 1959. XI + 182 und 1598. The Indian 
Society of Theoretical and Applied Mechanics. Kha- 
ragpur, India. 


M. Barner, Differential- und Integralrech- 
nung, Band I: Grenzwertbegriff, Differentialrech- 
nung. (Sammlung Göschen, Band 86/86a.) 1768. 
Berlin 1961. Walter de Gruyter & Co., Preis brosch. 
DM 5,80. 

W. Prager, Einführung in die Kontinuums- 
mechanik. 22858. m. 26 Abb. Basel und Stuttgart 
1961. Birkhäuser Verlag. Preis geb. DM 32,50. 


D. Meksyn, New Methods in Laminar Boun- 
dary-Layer Theory. XXIV + 294 S. Oxford/ 
London/New York/Paris 1961. Pergamon Press. 
Preis geb. 70 s. 


6. Hutarew, Regelungstechnik. Kurze Einfüh- 
rung am Beispiel der Drehzahlregelung von Wasser- 
turbinen. Zweite neubearbeitete Aufl. XII + 180 8. 
m.196 Abb. Berlin/Göttingen/Heidelberg 1961. Sprin- 
ger-Verlag. Preis geb. DM 21,—. 


H. Blasius, Lagerreibung. Hydrodynamische 
Theorie nach Sommerfeld und Michell, eine Darstel- 
lung der Grundlagen. 348. Hamburg 1961. Boysen 
& Maasch Verlag. .Preis brosch. DM 4,50. 

1.1. Goldman and V.D. Krivehenkov, Problems 
in Quantum Mechanics. 275 S. m. 33 Abb. Ox- 
ford/London/New York/Paris 1961. Pergamon Press. 
Preis geb. 50 s. 


NACHRICHTEN 


3. Internationaler Kongreß für Kybernetik 
(11.—15. 9. 1961 in Namen (Namur), Belgien) 


Ergänzend zur Vorankündigung auf S. 384 des vori- 
gen Jahrgangs wird auf Wunsch der Tagungsleitung 
folgendes mitgeteilt: 


Vortragsanmeldungen werden an die unten ange- 
führte Anschrift erbeten; für Vortragende ist die Teil- 
nahme gebührenfrei. Die Teilnahmegebühr beträgt für 
Mitglieder der AIC B.F. 200,—, sonst B.F. 400, —. 


Nähere Auskünfte erteilt: 


Association Internationale deCybernötique,a.s.b.1., 
Secrötariat, 13, Rue Basse-Marcelle, 
Namur (Belgique). 


DMYV-Tagung 1961 


Die diesjährige Tagung der Deutschen Mathemati- 
ker-Vereinigung wird vom 17. bis 24. September in 
Halle (Saale) durchgeführt. Im Programm ist u.a. 
eine Sektion über Wahrscheinlichkeitsrechnung und 
Statistik vorgesehen mit Vorträgen von B. W. GNE- 
DENKO, K. JAcoBs, A. N. KoLmocorovy und K. 
KRICKEBERG. Die Mitglieder der GAMM sind zu dieser 
Tagung herzlich eingeladen. Interessenten wollen sich 
bitte an das 

I. Mathematische Institut 
Halle (Saale) 
Universitätsplatz 8/9 


wenden, damit ihnen die Tagungsunterlagen zuge- 
schickt werden können 
Prof. Dr. O.-H. KELLER 
Vorsitzender der DMV 
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IFIP-Kongreß 1962 


Aufruf zur Anmeldung von Vorträgen 

Die IFIPS (International Federation of Infor- 
mation Processing Societies) wird vom 27. August bis 
zum 1. September 1962 in München einen Kongreß 
abhalten. $ 

Gegenstand des Kongresses werden alle Gesichts- 
punkte der Informationsverarbeitung und der digi- 
talen Rechenanlagen sein einschließlich folgender 
Teilgebiete: 


1. Informationsverarbeitung in der Wirtschaft („Bu- 
siness Information Processing“) 

Z.B. Datenverarbeitung in Handel, Industrie und 
Verwaltung (data processing in commerce, industry, 
and administration). 


2. Informationsverarbeitung in der Wissenschaft 
(„Seientifie Information Processing‘) 


Z. B. numerische Mathematik (numerical analysis); 
Rechenverfahren der angewandten Mathematik, Sta- 
tistik und Technik (caleulations in applied mathe- 
matics, statistics, and engineering); Datenreduktion 
(data reduction); Probleme der Unternehmensfor- 
schung (operations research). 


3. Zeitechte Informationsverarbeitung („Real Time 
Information Processing“) 


Z.B. Reservierungssysteme (reservation systems); 
Regelung durch Rechenanlagen (computer control); 
Verkehrsregelung (traffic control); Analog-Digital- 
Konvertierung (analog-digital conversion). 


4. Informationsspeicherung und -wiederauffindung 
(„Storage and Retrieval of Information“) 


Z.B. Speichervorrichtungen (memory devices); 
Bibliotheks-Kataloge (library catalogues). 


5. Sprachübersetzung und Sprachanalyse („Language 
Translation and Linguistie Analysis“) 


6. Digitale Nachrichtenübermittlung (‚Digital Com- 
munication‘) 


2. B. Verschlüsselung (encoding); Entschlüsselung 
(decoding) ; fehlerentdeckende und fehlerkorrigierende 
Codes für digitale Datenübertragung (error detecting 
and error correcting codes for digital data trans- 
mission). 

7. Künstliche Wahrnehmung und Intelligenz (,„Artifi- 
ceial Perception and Intelligence“) 

2. B. Erkennen von Formen (pattern recognition); 
biologische Modelle (biological models); lernende Ma- 
schinen (machine learning); Automatentheorie (auto- 
mata theory). 


8. Neuere Fortschritte bei Rechenanlagen („Advanced 
Computer Techniques“) 

Z. B. logischer Entwurf (logical design); logische 
Elemente (logical elements); Speichervorrichtungen 


(storage devices); überschnelle Rechenmaschinen 
(ultra high-speed computers); Programmierungs- 


methoden (program techniques); ALGOL. 
9. Ausbildung 


2. B. Auswahl und Schulung von Rechenmaschi- 
nen-Spezialisten (selection and training of computer 
specialists); Einführung von Nicht-Spezialisten in den 
Gebrauch von Rechenmaschinen (training of non- 
specialists in the use of computers); Informations- 
verarbeitung im Hochschullehrplan (information pro- 
cessing as a University subject). 
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10. Verschiedenes FE NOIR 

Z. B. die Ausbreitung des Gebiets der 
verarbeitung. int ee 

In jeder Kategorie sollen, soweit möglich, 
wendung digitaler Rechen 
Systementwurf, logisc 
Einzelteile behandelt werden. 


Er ws 


An 


Vortragsanmeldungen, begleitet von einer Zusam- er 
menfassung in englischer Sprache im Umfang von 500 


bis 1000 Worten, sollen bis . 
15. September 61° 


an den DARA-Programmausschuß für IFIP-Kon- 
greß 62 (Obmann: Prof. Dr. F. L. Bauer, Institut für 
Angewandte Mathematik der Johannes-Gutenberg- 
Universität, Mainz, Jacob-Welder-Weg 7; stellver- 
tretender Obmann: Dr. H. Bırrıng, Max-Planck- 
Institut für Physik und Astrophysik, Institut für 
Astrophysik, Abt. Numerische Rechenmaschinen, 
München 23, Aumeisterstr.) eingereicht werden. 


Die Vortragsanmeldungen werden auf Grund der 
Zusammenfassungen vom internationalen Programm- 
Komitee der IFIPS geprüft werden. Die Autoren der 
angenommenen Zusammenfassungen werden bis zum 
März 1962 eine Aufforderung erhalten, ihre vollstän- 
digen Manuskripte (in Französisch oder Englisch) dem 
Programm-Komitee zur endgültigen Prüfung vorzu- 
legen. 

Zusätzlich zu den ausgewählten Vorträgen werden 
Vorträge auf Einladung, Symposien und vorbereitete 
Diskussionen (‚‚panel discussions‘‘) stattfinden. An- 
regungen hierfür mögen ebenfalls an den DARA- 
Programmausschuß gerichtet werden. 


BERICHTIGUNGEN 


Zu H.Störmer, Über einen Exponentialverteilungs- 
test, ZAMM 40 (1960), 8. T 6—T 97: 


Im oben zitierten Vortragsauszug muß es heißen 
2v—1 
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A= [mas 


Fe 1 
2(n —1)| ' 
Mit dem ursprünglich angegebenen Bereich A erhält 
man eine größere statistische Sicherheit als S. 
H. STÖRMER 


Zu J.Wasowski, Zwei Lehrsätze in der Operatoren- 
rechnung, ZAMM 41 (1961), S. 9—17: 

Nach Mitteilumng des Verfassers sind folgende 
Korrekturen nachzutragen: 

l. In der Zeile unter Formel (1.4a) ist der Hinweis 
auf die Fußnote durch den Hinweis auf die Literatur- 
stellen [1], [2], [4] zu ersetzen. 

2. Formel (3.1) lautet: 


u er 
e)= 3 —pP0 +). 
p=0 p: 


D 


3. Auf S. 16 lautet der erste Summand in der vier- 
ten Zeile von unten: 6 a9). 


Herausgeber und Hauptschriftleiter: Prof. Dr.-Ing. habil, H. Heinrich. D 27 ied 
geber und tsch eiter: Prof. Dr.-Ing. Re ‚ Dresden A 27, Friedrich-Hegel-Str. 31. Verlag: Akademie-Verla. 
GmDE, her W 8, Leipziger Straße 3—4; Fernsprecher: 220441. Telex-Nr. 011773. Postscheckkonto: Berlin 35021. A ee nee 
eftes: 1009/41/6. Die Zeitschrift für Angewandte Mathematik und Mechanik erscheint monatlich. Bezugspreis: vierteljährlich DM 15,—. 


Zuzüglich Bestellgeld. Doppelheft DM 10,—, 


Abbestellungen können nur bis 4 Wochen vor Quartalsende anerkannt werden, andernfalls 


wird das folgende Quartal noch geliefert, Veröffentlicht unter Lizenznummer ZLN 5011 des Ministeriums für Kultur. Gesamtherstellung: 
VEB Druckerei „Thomas Müntzer‘‘ Bad Langensalza (V/12/6) (1). Printed in Germany. 


Diesem Heft liegen zwei Prospekte des Springer-Verlages Berlin bei. 
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henanlagen, ung, 
her Entwurf, Ausrüstung und 


EREL WINKLER 


Elektronische Analogicanlagen 


(Bickteonisches Rechnen und Regeln, Band 2) 


1961. X. 242 Seiten — 172 Abbildungen — gr. 8° — Lederin DM 29,— 


In dem Band geht der Verfasser auf das Wesen und die Bedeutung von Bes und Modellen 
ein und erläutert die Problemerfassung mittels Wirkschaltbilder, Blockschaltbilder, Struktur- . 
bilder und Koppelpläne, um letzthin analogieelektronische Lösungswege beschreiten zu kön- 
nen. Er gibt eine Einteilung der elektrischen Analogieanlagen, eine Zusammenstellung aller 
linearen und nichtlinearen Analogieelemente und beschreibt ihr Zusammenwirken. Ferner be- 
handelt er systematisch und ausführlich die mannigfachen Anwendungen elektronischer 
Analogieanlagen mit dem Ziel, den Leser zum Einsatz analogieelektronischer Methoden auf 
seinem Spezialgebiet anzuregen, um den Fortschritt zu beschleunigen, Zeit und kostspielige 
Versuchsaufbauten zu sparen. 

In dem zweiten Teil des Buches schildert der Autor die physikalische und technische Seite, 
die Entwicklungs- und Konstruktionsmerkmale der elektronischen Analogieanlagen, ihre 
Analogieelemente und Baugruppen: Operationsverstärker und elektrische ‚Netzwerke als 
lineare Analogieelemente, nichtlineare Analogieelemente und Funktionsempfänger, Hilfs- und 


Zusatzeinrichtungen wie Meßwandler für den Simulationsbetrieb und schließlich Typen- 
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beispiele kompletter elektronischer Analogieanlagen. Abschließend vermittelt er einen histori- 
schen Überblick. Besonders hervorzuheben ist der im Anhang gegebene sehr umfangreiche 
Literaturnachweis, der das einführende Werk wertvoll.ergänzt. Das Buch ist das erste und 
bisher einzige auf dem Gebiet der Analogieelektronik, das von einem deutschen Verfasser er- 


scheint. 


Bestellungen durch eine Buchhandlung erbeten 


BER LTIN 


Karl- ER Leipzig am Mi t 
gehalten hat. Erstmalig in der deutschen] 
iR auf der Grundlage des d 


wobei vor allem die ntyicktunglinicn heran e 


Beziehungen nachekor ieh. 


Bestellungen nimmt jede Buchhandlung entgegen 


EB 


B.G.Teubner Verlagsgesellschaft Leipzi 


